Properties of neutral mesons in a hot and magnetized quark matter 
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The properties of noninteracting a and ir° mesons are studied at finite temperature, chemical po- 
tential and in the presence of a constant magnetic field. To do this, the energy dispersion relations 
of these particles, including nontrivial form factors, are derived using a derivative expansion of the 
effective action of a two-flavor, hot and magnetized Nambu-Jona-Lasinio (NJL) model up to second 
order. The temperature dependence of the pole and screening masses as well as the directional 
refraction indices of magnetized neutral mesons are explored for fixed magnetic fields and chemical 
potentials. It is shown that, because of the explicit breaking of the Lorentz invariance by the mag- 
netic field, the refraction index and the screening mass of neutral mesons exhibit a certain anisotropy 
in the transverse and longitudinal directions with respect to the direction of the external magnetic 
field. In contrast to their longitudinal refraction indices, the transverse indices of the neutral mesons 
are larger than unity. 
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I. INTRODUCTION 

The study of the states of quark matter under ex- 
treme conditions has attracted much attention over 
the past few years. Extreme conditions include high 
temperatures and finite baryonic chemical potentials 
as well as strong magnetic fields. The latter is re- 
sponsible for many interesting effects on the prop- 
erties of quark matter. Some of the most impor- 
tant ones are magnetic catalysis of dynamical chiral 
symmetry breaking dHl], that leads to a modifica- 
tion of the nature of electroweak [J], chiral and color- 
superconducting phase transitions ja-llOl]. production 
of chiral density waves [ll| , chiral magnetic effect [l^ , 
and last but not least inducing electromagnetic super- 
conductivity and superfluidity In this paper, we 
will focus on the effect of constant magnetic fields on 
the properties of neutral and noninteracting mesons in 
a hot and dense quark matter. In particular, the tem- 
perature dependence of meson masses as well as their 
direction-dependent refraction indices 1 and screening 
masses will be explored in the presence of various fixed 
magnetic fields. The largest observed magnetic field 
in nature is about 10 12 — 10 13 Gaufi in pulsars and 
up to 10 14 — 10 15 Gaufi on the surface of some mag- 
netars, where the inner field is estimated to be of or- 
der 10 18 - 10 20 Gaufi There are also evidences 
for the creation of very strong and short-living mag- 
netic fields in the early stages of non-central heavy 
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1 The term "refraction index" is used in [14( for pions modified 
by the matter (qua sipions). Although, the same terminology 
is also used in [\M , the definitions of refraction index in |14ll 
and |15|I are slightly different, as will be explained later. 



ion collisions at RHIC [17|, [l8| ■ Depending on the col- 
lision energies and impact parameters, the magnetic 
fields produced at RHIC and LHC are estimated to 
be in the order eB ~ 1.5 to 2 , corresponding to 0.03 
GeV 2 for m n = 138 MeV, and eB ~ 15 to 2 , cor- 
responding to 0.3 GeV 2 , respectively [lj|. 2 On the 
other hand, it is known that the quark-gluon plasma, 
produced in high-energy heavy-ion collisions, passes 
over many stages during its evolution. The last of 
which consists of a large amount of hadrons, includ- 
ing pions, until a final freeze-out Thus, the pres- 
ence of a background magnetic field created in heavy 
ion experiments may affect the properties of "charged 
quarks" in the earliest stage of the collision and al- 
though the created strong magnetic field is extremely 
short living and decays very fast [HI, [l9| , it may affect 
the properties of the hadrons made of these "magne- 
tized" quarks. Even the properties neutral mesons 
may be affected by the external magnetic field pro- 
duced in the earliest phase of heavy-ion collisions. In 
the present paper, we do not intend to go through the 
phenomenology of heavy-ion collisions. Our computa- 
tion is only a theoretical attempt to study the effect 
of external magnetic fields on "magnetized" neutral 
mesons, that, because of the lack of electric charge, 
do not interact directly with the external magnetic 
field. Our study is indeed in contrast with the recent 
studies in 2(J 21 1, where chiral perturbation theory 
is used to study the effect of external magnetic fields 
on the pole and screening masses as well as the decay 
rates of charged pions interacting directly with the 
external magnetic field. 

There are several attempts to study the effect of 
temperature and chemical potential on the properties 



Note that eB = 1 GeV 2 corresponds to B ~ 1.7 X 10 20 Gaufi. 



pions in a hot and dense medium, in the absence of 
external magnetic fields EE H2411. In Q, the 
energy dispersion relation of the so-called "quasipi- 
ons" (or pions modified by the matter) is introduced 
by 



J (P) 



u 2 p 2 + m\. 



(LI) 



Here, u(T) is the temperature-dependent refraction 
index (also called "mean quasipion velocity" EH), and 
is the pole mass of the pions. To determine m^, 
one can either start from the Lagrangian density of a 
linear er-modcl including four-pion interaction or use 
the chiral perturbation theory Lagrangian within cer- 
tain approximation. Considering the pion (one-loop) 
self-energy of the model, and computing, in particu- 
lar, its pole, it is possible to determine the pion pole 
mass (at one- loop level). As concerns the screening 
mass of pions, m*, it is related to through the 
relation m£ = m^/v^, where v n is the pion velocity 
|22| . As it is shown in 22 1, the velocity of massless 



pions is in general given by 



2 2 2 _ -R- e /l 2 



(1.2) 



where w = po is the energy, p = |p| is the absolute 
value of pion three momentum, and f* and /* are tem- 
poral and spatial pion decay constants, respectively. 
As it turns out, at zero temperature, because of rel- 
ativistic invariance, /* = and therefore v„ = 1. 
At finite temperature, however, since a privileged rest 
frame is provided by the medium, relativistic invari- 
ance does not apply anymore, and, as it is shown in 
[lU, "cool" pions propagate at a velocity v„ < 1. 
Moreover, it is shown in [22j |. that for approximate 
chiral symmetry, the Gcll-Mann, Ookes and Renner 
(GOR) relation between the pion mass and the 
pion decay constant /„. still holds at finite tempera- 
ture, except that instead of f v , the real part of /* 



enters the GOR relation, i.e. 



2m o ('0-i/j) 

(He/* 



Let us 



also notice that at finite temperature and in the ab- 
sence of external magnetic fields, no distinction is to 
be made between neutral and charged pion masses. 

Nontrivial energ y d ispersion relation of mesons is 
also introduced in jl5| | and (24|. In EH, using the 
definition of the group velocity, a momentum de- 
pendent "refraction index" n(p) is defined for pions 
by the ratio of the group velocity in matter and in 
vacuum, n(p) = v^/v gr . Here, the matter pion 

group velocity is defined by v gr = with po = 

[n^iT,^ 2 + M 2 (T, /i)] 1 / 2 , and the vacuum pion 
group velocity is defined by w™° = The mo- 

mentum dependent refraction index is therefore given 
by n(p) = (^p£k^J n. It is argued that since for fi- 
nite temperature T and chemical potential /i, we al- 
ways have both n > 1 and 4j§c > 1 for all values 
of p, the index of refraction developed by the pion 



medium at finite T and fi is always larger than unity 
EH . Let us notice that the definition of the refraction 
index n in J15I is slightly different from what is used 
in 14 [: In |15| . n _1 appearing in the dispersion rela- 
tion po = [n _1 (T,/i)p 2 + M 2 (T, [i)] 1 / 2 is the same as 
u 2 appearing in the dispersion relation (11.11) from EH • 
In the latter, u = n^ 1 / 2 is called refraction index. 3 
Having this in mind, it turns out that the results pre- 
sented in El. coincides with those obtained in [24| . 
Here, the quantity u appears as in (l4T |. in the pion 



energy dispersion relation, 



-(p 2 + m 2 ), and is 



termed "velocity" , although the authors mention that 
u is the pion velocity only when m = 0. Here, m is the 
screening mass. The pion pole mass is then defined by 
m p — um. Using scaling and universality arguments, 
the authors predict that "when critical temperature 
is approached from below, the pole mass of the pion 
drops despite the growth of the pion screening mass. 
This fact is attributed to the decrease of the pion ve- 
locity near the phase transition" [24j. 

As concerns the effect of external magnetic fields 
on the low energy properties of QCD, in [25], the GOR 
relation between the neutral pion mass m w o and its 
decay constant / T o, is shown to be valid in the first 
order a chiral perturbation theory in the presence of 
constant and weak magnetic fields, whose Lagrangian 
includes, in particular, (7? 2 ) 2 self-interaction terms. 
This method is also used recently in [IE Hl| to deter- 
mine the pion thermal mass and the pion decay con- 
stants in the presence of a constant magnetic field and 
at finite temperature. It is shown, that the magnetic 
field gives rise to a splitting between m^o and rn 7r ± as 
well as /^o and f n ± . The pion decay constants f n o and 
f v ± are computed by evaluating the matrix elements 
(0|A°|7T°) and (0|Aj|7r T ), respectively. However, no 
distinction between the temporal (fi = 0) and spatial 
(/i = 1, 2, 3) directions is made. 

In the present paper, we will mainly focus on non- 
trivial energy dispersion relations of noninteracting a 
and 7? mesons, arising from an appropriate evaluation 
of the one-loop effective action of a two-flavor NJL 
model in a derivative expansion up to second order. 
Our method is therefore different from the method 
used in [IE Hl|, and involves, in contrast to [IE [2lj . 
the effect of external magnetic fields on charged quarks 
from which the mesons are built. This will give us the 
possibility to explore the effect of external magnetic 
fields on neutral mesons at finite temperature and 
chemical potential. Using the method originally intro- 
duced in (E [26| for a single flavor NJL model, we will 
arrive at the effective action of a and 7? = (tti, 7T2, 7^) 



3 In the present paper, we have adopted the terminology used 
in 0. 



2 



mesons, 



r eff [cr,7?] = r off [cr ] 



1 



-- J d d xa(x) (Ml + G^dl)a{x) 

3 

~E/ ^ **0"0 (M| + .7^)^(z), (1.3) 

including nontrivial meson squared mass matrices 
{Ml, Ml) and form factors {O^ ,T» V \ and leading 
to the energy dispersion relations of a and tt mesons 



2 i 2 

1 m„.. 



El = J>»ft) 



(1.4) 



Here, the pole masses (m 2 , m|) and refraction indices 
(u CT , Ujf) of the mesons are defined by 



Re[M 2 ] 
Re[£ 00 ] 



Re[M|J 
Re[(-F° )* 



and 



/ Re[g"] 
v Re[£ 00 ] 



1/2 



1Tf 



Ue[(^ 00 )< 



1/2 



for all space directions i — 1,2,3 and isospin indices 
£ = 1,2,3. These quantities can be computed using 
the one- loop effective action of a two- flavor N JL model 
at finite temperature T, chemical potential fi and for 
a constant magnetic field B 7 according to the formal- 
ism presented in 0, [2rt. Using the definition of the 
screening mass from [24|, the screening masses of a 

and 7? mesons, 



and m9 are given by 



7M'- 



and 



(*) : 



V« = 1,2,3, 



respectively. Later, we will, in particular, show that in 
the presence of a uniform magnetic field, directed in a 
specific direction, the refraction indices and screening 
masses in the transverse and longitudinal directions 
with respect to the direction of the background mag- 
netic field will be different. 

The organization of this paper is as follows. In 
Sec. mi we will generalize the method introduced in 
0] to a multi-flavor system, and will derive the ef- 
fective action (11.31) . using an appropriate derivative 
expansion up to second order. In Sec. IIII1 we will 
determine the one-loop effective potential of a two- 
flavor NJL model including (a, 7?) mesons. In Sec. 
TV\ the squared mass matrices (M 2 ,Af| ) and ki- 
netic coefficients (G^, F^ v ) corresponding to neutral 
mesons a and 7r° will be analytically computed at 
finite (T, zx, eB) and up to an integration over p%- 
momentum as well as a summation over Landau lev- 
els. In Sec. IV Al we will first use the one-loop effec- 
tive potential, evaluated in Sec. IIII1 to explore the 



phase portrait of the model. Here, the effect of mag- 
netic catalysis [l|, and inverse magnetic catalysis 
d, [I?} on the critical (T, /i, eB) will be scrutinized. 
Performing numerically the remaining ^-integration 
and the summation over Landau levels from Sec. IIV1 
we will present, in Sec. IV Bl the T-dependence of 
(M 2 ,M| ) and (G^ , J 7 **") for various fixed magnetic 
fields and at fi = 0. Using these results, the In- 
dependence of pole masses of neutral mesons as well 
as their directional refraction indices and screening 
masses will be determined in Sec. IV CI for various 
fixed eB = 0.03,0.2,0.3 GeV 2 . We will in particular 
show that, for non-vanishing magnetic fields, the re- 
fraction index of noninteracting mesons in the longitu- 
dinal direction is equal to unity, while their transverse 
refraction index is larger than unity. Let us notice 
that since the mesons are massive, this does not mean 
that magnetized mesons propagate with speed larger 
than the speed of light. 4 The observed anisotropy 
in the meson refraction indices is because of the ex- 
plicit breaking of Lorentz invariance by uniform mag- 
netic fields. The same anisotropy is also reflected in 
the screening masses of neutral mesons in the longi- 
tudinal and transverse directions with respect to the 
direction of the background magnetic field. We will 
plot the T-dependence of mesons screening masses for 
various fixed eB and /i, and will show that, in the 
transverse directions, they are always smaller than 
the screening masses in the longitudinal direction. 
Motivated by recent experimental activities at RHIC 
and LHC, we will only consider the effects of rela- 
tively weak and intermediate magnetic field strength 
(eB = 0.03,0.2,0.3 GeV 2 ). As concerns the effect of 
stronger magnetic fields, we will show that they lead 
to certain instabilities at low temperature. Our re- 
sults for eB — 0.5,0.7 GeV 2 are consistent with the 
main conclusions presented recently in [29], where a 
single flavor NJL model is studied in 2 + 1 dimensions 
in the presence of a strong magnetic field and at fi- 
nite temperature. A summary of our results will be 
presented in Section IVT1 



II. 



MATHEMATICAL TOOL: DERIVATIVE 
EXPANSION OF THE QUANTUM 
EFFECTIVE ACTION 



Let us consider a theory containing real scalar 
fields (ipo, <fi" ' i <Pn-i) = $j whose dynamics are de- 
scribed by the effective action r otf [$]. Using an appro- 
priate derivative expansion, and, inparticular, gener- 
alizing the method introduced in 0, [2f| to a multi- 
flavor system, we will derive, in this section, the en- 
ergy dispersion relations of tpe, £ = 0, • • • , iV — 1. Us- 



4 The effect of constant magnetic fields on the propagation of 
massless particles is recently discussed in |2S| . 



3 



ing the energy dispersion relation, the pole and screen- 
ing mass as well as the directional refraction index 
corresponding to tpi, I = 0, • • • , N — 1 will be defined. 

Let us start by expanding around an x- 

independent configuration $o 5 



$(x) = $ + $(x). 



(11.1) 



Plugging pi.ip in the effective action, we arrive first 
at 



r oft [$] = r cff [$ ] + / 



d x- 



5<pi(x) 



tpi(x) 



-i [d d xd d y 

2 J 6ipi(x)5ipj(y) 



ipi(x)ifj(y) + 



(11.2) 



Assuming that $o describes a configuration that min- 
imizes the effective action, the second term in (|II.2[) 
vanishes. Using then the Taylor expansion 



9(y) = §(x) + z"d^(x) + -z»z v d^d v ®(x) + 



(II.3) 



with z = y — x, and neglecting the terms linear in z, 
we get 

r eff [$] = r cff [$ ] - \J ddx M%[^o]ipi{x)(pj{x) 



+ ^ I ddx X%l®o}d^(x)d^(x) + • • • , (II.4) 



where the summation over i, j = 0, • • ■ N— 1 is skipped. 
In (|II.4[) . the "squared mass matrix" M 2 ^ and the "ki- 



netic matrix" Xu are given by 



M?-[$ ] = - J d d i 



s 2 r of{ 



6tpi(0)Sifj(z) 



'!>„ 



(II.5) 
.(II.6) 



The above derivative expansion of r off [$] from (|II.4j) 
can alternatively be given as 



r otf [$] = fd d x(-vm 



const., that spontaneously breaks the O(N) symmetry 
of the original action. Using (111.51) . or equivalcntly 



M 2 im [$ ] 



d d z 
d d z 



s 2 r c!! 



8<p (z)8ip (0) 

<5 2 r cff 



6tpi(z)8tp m (0) 



(II. 



V i,m > 1, it is possible to determine the squared 
mass matrices corresponding to the collective modes 
ipo,ipi,--- , <£jv-i. To determine the kinetic part of 
the effective action, we use, as in [2j, the Ansatz 



(11.9) 



$2 ' 

V i,j = Q,l,>--N- 1. Here, <S> 2 = ^? and 

xf/[$o] = xf/[$o], appearing in pL4]). Plugging 
(111.91) in pi.7p . the kinetic part of the effective La- 
grangian density including two derivatives is given by 

1 „ 

£fc = ^(F^ v )i j d^ l d lJ f ] + -p- (cpid^ipi) ((pjd u cpj) . 

(11.10) 

To determine the form factors F± v and i^", or at 
least a combination of these two form factors, we will 
use the definition of r* ff = J d d xCu , as a part of the 
effective action including only two derivatives • We 
get 



2pfc 



s 2 r 



Sip o (x)S(fo(0) 
with Q^ v = [(F? v ) 



d fl dj d (x), (11.11) 



2F% V \ , and 



2 r k 



S 2 T 



Sip e (x)Sip m (0) 



d^8 d (x), (11.12) 



V £, m > 1, where (T^)tm = | [(F^W + (Ff)™/] 
From (|ILTT|) and pTT2)| we have 



SH^o] = -r / d d zz»z» - ; ff 



(0) 



6<pe(z)8tp m (Q) 



*0 

(11.13) 



tfiW^fc) 8 ^ (z) + • • • ) , (II.7) V £, to > 1. Comparing the above relations with x 

from (|II.6|) . it turns out that Xoo = and Xtm 



where, all non-derivative terms in (|II.4j) are summed 
up into the potential part of the effective action U[$], 
and the terms with two derivatives yield the kinetic 
part of the effective action, proportional to Xij- To 
have a connection to the example that will be worked 
out in the subsequent sections, let us assume a fixed 
configuration for $o = (<A)(o) > 0> 0> ' " ' > 0) , with <po(o) = 



tin 



(F^)i m y £,m > 1. Assuming then (M 2 ) 
-(M 2 )mt, and (T^) tm = -(T^) mU V l± to and 



, to > 1, and denoting A^qq by Mg, as well as M. 



5 This will be shown in our specific example in the subsequent 
sections. 



4 



by M\ for £ = 1, • • • ,N—1, the effective action (|IL4|) 
simplifies as 

r cff [<f] = r off [$ ] - \ J d d xt p (m 2 + g^dl) $ 

N-l 

— J2 / d*x<p t [M2 + {F^) u dl)w. (11.14) 

Here, we have used the fact that Q^ v and J 7 ^ are di- 
agonal, i.e. Q» v = QWg^ as well as = jw^. 
The same relations are shown to be valid in a single- 
flavor case [2|. From (jll. 14[) . the general expressions 
for the energy dispersion relation of noninteracting 
<pt,£ = 0, ■ ■ • ,N — 1 fields can be determined. For 
£ — 0, we have 

Kc = (^V? + S 22 ^ + ^ 33 Pi + M 2 ) , (11.15) 
and for W > 1, we have 



1 



• [(J-")„ V \ + (^ 22 )« p£ + (^ 33 )« P 2 + M 2 ] . 

(11.16) 



Using the above energy dispersion relations, the pole 
masses of free tpo and if£, I > 1 are given by 



l2 _ M 



and rri£ 



M 2 



-, (11.17) 



respectively. The screening masses Trip, and "direc- 
tional" refraction indices up of noninteracting ip^, t — 
0, 1, • • • , N — 1 fields in the i-th directions (i = 1, 2, 3) 
are defined by 



(i) "10 



^, where (n«) 2 = ^, (11.18) 



.(0 

for i — 0, as well as 

(i) _ 



where (it], ) 



KK2 = P 7 " 
£ 7 (J-oo) 



for £ > 1 [see Sec. |V]for more details on the definition 
of screening masses and refraction indices] . 

In the present paper, we will use the above dis- 
persion relations, to describe the properties of non- 
interacting a and 7? mesons in a hot and magnetized 
medium. We will focus, in particular, on a and 7T3 
mesons. The latter will be identified with the neutral 
pion, 7T3 = 7T°. To do this, we will first consider, in 
the next section, a two-flavor NJL model including ap- 
propriate four-fermion interactions. Defining the me- 
son fields a and tt in terms of fermionic fields, and 
eventually integrating the fermions in the presence 
of a constant magnetic field, we arrive at the one- 
loop effective action r eff [(T, 7?], describing the dynam- 
ics of magnetized meson fields. We will spontaneously 



break the chiral symmetry of the original theory, by 
choosing a fixed configuration (<to,7?o) = ( cons t-jO), 
that minimizes r eff [<7, 7?]. Using then the formalism 
described in the present section for the specific case 
of N = 4, and identifying (po with the a-meson and 
tpe, £=1,2,3 with the pions 717, £ = 1,2, 3, we will de- 
termine the temperature dependence of the pole and 
screening mass, as well as the directional refraction 
indices of noninteracting neutral a and ir° mesons at 
finite temperature and in the presence of a constant 
magnetic field. We will postpone the discussion on 
the properties of charged and magnetized pions to a 
future publication [30j . 



III. ONE-LOOP EFFECTIVE POTENTIAL 
OF A TWO-FLAVOR NJL MODEL AT FINITE 

(T,/x,eB) 

In this section, we will determine the one-loop effec- 
tive potential corresponding to a two-flavor magne- 
tized NJL model at finite temperature and density. 
The minima of this effective potential will then be 
used in the subsequent sections to determine the ki- 
netic coefficients and mass matrices corresponding to 
neutral a and 7r° mesons. 

Let us start by introducing the Lagrangian density 
of a two-flavor gauged NJL model in the presence of 
a constant magnetic field 

C = $(x)(i^D fl -m )i>(x) + G {$(x)iP(x)} 2 



+$(x)i l5 TiP(x)} 2 } - -F^F^. 



(III.l) 



Here, the fermionic fields ip'j carry apart from the 
Dirac index, a flavor index / S (1,2) = (u,d) and 
a color index c 6 (1,2,3) = (r,g,b). In the chiral 
limit m — > 0, this implies the SUl(2) x SUr(2) chi- 
ral and 5/7(3) color symmetry of the theory. The 
isospin symmetry of the theory is guaranteed by set- 
ting m u = m,4 = m,Q. The covariant derivative in 
-, (11.19) (|HL1~J) is defined by D» = <9 M + ieQAf*- , where Q = 
diag (2/3, — 1/3) is the fermionic charge matrix cou- 
pled to the U(l) gauge field Aff*', and f = (n, r 2 , r 3 ) 
are the Pauli matrices. Choosing, the vector poten- 
tial A ^ L in the Landau gauge A™*- = (0,0,^,0), 
pil.ip describes a two-flavor NJL model in the pres- 
ence of a uniform magnetic field B = Be$, aligned in 
the third direction. The field strength tensor F^ u is 
defined as usual by F^ u = d^A^-, with A^f*- fixed 
as above. As it turns out, the above Lagrangian is 
equivalent with the semi-bosonized Lagrangian 



C sb = ip(x) (ij^D^ - m )ip(x) - (a + i~/ 5 T ■ Tr)ip 
(a 2 + tt 2 ) B 2 
AG 2 ' 



(III.2) 



where the Euler-Lagrange equations of motion for the 



5 



auxiliary fields lead to the constraints 

a{x) = -2G4>(x)ip(x), 
7?(x) = —2Gip(x)ij5Tip(x). 



(111.3) 



To determine the one-loop effective action correspond- 
ing to (jlll.f [) as a functional of a and tt, the fermionic 
fields ip and ?/> in (|III.2[) are to be integrated out. Us- 



ing 



ir eff [<r,7?] 



VipTty exp [i / d 4 x C sb , (111.4) 



the one-loop effective action T aB is then given by 

r ea [a, tT] = r£ } [a, tt] + r$ [a, tt] , (111.5) 

where the tree level part, , and the one-loop part, 
, are given by 



rS ) [a, 7 rl = - / d 4 x 



and 



riff [tr, tt] = -iTr {cfsx} ]n[iS Q l (<T,i?)]. (HI.7) 



Here, to = too + cr(a;) and 

iS Q \a, tt) = i^Dp - (to + i 7 5 f • 7?) , (111.8) 

is the inverse fermion propagator. To determine 
rlf f ^ [tr, 7r] , let us assume a constant and fixed config- 
uration (tTo,7?o) = (const., 0) for the collective modes 
(tr, tt), that breaks the SUl{2) x SUr(2) chiral sym- 
metry of the original action in the chiral limit. Only 
in this case, to can be replaced by the constant con- 
stituent quark mass to = too + tro, where (To =const. 
The one-loop effective potential is given by evaluating 
the trace operation in pil.7|) . that includes a trace 
over color c, flavor /, and spinor s degrees of freedom, 
as well as a trace over a four-dimensional space-time 
coordinate x. Following the standard method intro- 
duced e.g. in 0, and after a straightforward computa- 
tion, the one-loop part of the effective action rl f y [tro] 
reads 

rl g ) [t7 ] = -6i Yl lndet s [££-pg], (HI.9) 



where the energy of a charged fermion in a constant 
magnetic field is given by 



= J 2\qeB\p + pl + m 2 . (111.10) 



Here, the Ritus four-momentum 



p q = (p ,0, -Sgn(qeB)^/2\qeB\p,p 3 ), (III.ll) 



arises from the solutions of Dirac equation in the pres- 
ence of a constant magnetic field (see [U |32[ for 
more details on the Ritus Eingenfunction method). 
In 1 0[) . p labels the corresponding Landau lev- 
els appearing in the presence of a uniform mag- 
netic field. Performing the remaining determinant 
over the coordinate space in (lHL9|) leads to the ef- 
fective (thermodynamic) potential defined by 
filff = — V _1 rl f V, where the factor V denotes the 
four-dimensional space-time volume. The final form 
of olf f ^ is then determined in the momentum space, 
where the effect of finite temperature and chemical 
potential is introduced by replacing po in (|III.9j) with 
Po = iuj n — /i. Here, the Matsubara frequencies Lu n 
are defined by w„ = (2n + 1)itT. Using the standard 
replacement 



d A p ,, 



\qeB\ 



+00 +00 P+OC £ 

E E"? / o"T f( iu} n ~ M)P)P3), 

,_ ^ „- n J — CO 

(111.12) 



n— — 00 p— 



with p labeling the Landau levels and (3 = T _1 , and 
after summing over the Matsubara frequencies n, the 
(one-loop) effective potential of the model reads 



(i) 



-3 E 



\qeB\ 



+00 p+00 j 
x E «p / % {fM q + In (l + e -^+^>) 

+ ln(l + e-^ E "-^Y (111.13) 



Here, a p = 2 — 5 p0 is the spin degeneracy factor. As it 

turns out, the above expression for Cl[ B consists of a 
(T, /z)-independent and a (T, /^-dependent term. The 

(T, ^-independent part of Cl[g is divergent and is to 
be appropriately regulated. In the Appendix, we have 
followed the method presented in [33| , and shown that 
the (T, /x)-independent part of olff is given by (|A.12I) . 
Adding this part to the tree level part of the effective 
potential, (jIII.6l) . as well as to the (T, ^-dependent 

part of filff* , we arrive at the final expression for the 
one-loop effective potential of a two-flavor NJL model 
at finite (T, /z) and in the presence of a uniform mag- 
netic field aligned in the third direction 



G 



3 f „, / A + VA 2 + m 2 S 



4tt 2 



m In 



A(2A 2 + m 2 )VA 2 +m 2 



E 



\qeB\ 



ge{ |,-i} P =o 



+00 /-+00 J 

^ ap y glln^ + e-^+^+ln^ + e-^-"))}. (111.14) 



Here, x q = 2 \qeB \ ' ^ * s an a PP r0 P r i a te ultraviolet 
(UV) momentum cutoff, = ^ I 



and 



is given in (|III.10|) . In Sec. |V1 after fixing a num- 
ber of free parameters, such as the coupling G and the 
UV cutoff A, the global minima of f2 eff (m; T,fj,,eB) 
will be determined numerically. They will be then 
used to determine the squared mass matrices M 2 and 
M 2 and the form factors (kinetic coefficients) Q^" 
and (J-" Ml/ )33, corresponding to the neutral mesons a 
and 7T°, at finite (T,/j) and eB. 



IV. EFFECTIVE KINETIC PART OF THE 
ONE-LOOP EFFECTIVE ACTION OF A 
TWO-FLAVOR NJL MODEL AT FINITE 

(T, fi, eB) 

In the previous section, the one-loop effective poten- 
tial of a magnetized two-flavor NJL model at finite 
(T, /i.) is computed by evaluating the trace operation 
in (|III.7P for a fixed field configuration $ = (00 > tto) = 
(const., 0), which is supposed to minimize the one- 
loop effective potential (jlll. 14[) of the model. In the 
next two sections, we will compute the squared meson 
mass matrices and form factors of the effective kinetic 
part of the one-loop effective action corresponding to 
neutral mesons a and tt°. This computation includes 
an analytical and a numerical part. In this section, 
after reformulating the general derivation presented 
in Sec. UTJ and making it compatible with our case 
of magnetized two-flavor NJL model, we will present 
the analytical results of the squared mass matrices 
(M 2 ,M 2 ) and form factors {Q^ v ,T^ V ) for neutral 
mesons up to a one-dimensional integration over p 3 - 
momentum and a summation over Landau levels p. 
They shall be performed numerically. The results 
of the numerical computation will be presented in 
Sec. El where we explore the (T, /x, eB) dependence of 
(M 2 , M 2 ) and (Q^ , J 7 ^). Using these quantities the 
pole and screening masses of free neutral mesons and 
their directional refraction indices will be determined 
for various (T, jj,, eB). 

As we have described in Sec. HH our goal is to 
bring the effective action of a two-flavor NJL model 



including (a, 7?) mesons, in the form 

r cff [a, 7?] = r cff [cr ] 

' f d d xa{x){Ml+g^dl)a{x) 



"2 E / d " x (Ml+T^dl) u Tt t (x),(W.l) 

which is valid in a truncation of the derivative ex- 
pansion of the full effective action r Bft [<r, n] up to two 
derivatives. According to (|II.8j) . the squared mass ma- 
trices of neutral mesons, a and 7r°, are given by 6 



M, 



2 = - / d 4 - 



(M|)33 = - J 



Sa(0)Sa{z) 



(c-0,0) 



5n 3 (0)57r 3 {z) 



(IV.2) 



(0-0,0) 



and, according to (|II.13|) . the form factors of the effec- 
tive kinetic part of the effective action, corresponding 
to a and 7r°, read 



33 



2-nfc 



d 4 zz^z ] 
d i zz^z" 



6 2 T 



5a(0)da{z) 
S 2 Tt 



(0-0,0) 



57r 3 (0)6n 3 (z) 



. (IV.3) 

(o-o.O) 



To simplify our notations, we will denote in the 
rest of this paper, the mass squared matrix (Af2) 33 
from (|IV.2[) corresponding to ir° by M 2 . Similarly, 



X', 



will be denoted by T^ v . Whereas the mesons 



squared mass matrices at zero temperature and chem- 
ical potential are given by plugging the effective action 
fnT^ -fnO ]) in (1TV21) and read 



M 2 a 
M 2 



^ - i J d 4 ztr sfc [S Q (z, 0)Sq(0, z)} , (IV A) 



2G 



1 J d 4 ztr sfc [S Q (z,0)T 3l 5 S Q (0,z) 7 5 T 3 ], 

(IV.5) 



6 Here, the third component of 7? is identified with 7T , i.e. 
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the form factors (|IV.3|) arise by replacing r* f with the 
one-loop effective potential T^} from (|III.6P - (|III.8|) . 

= -1 J d 4 zz»z"tr sfc [S Q (z,0)S Q (0,z)], (IV.6) 

J*" = '-f d 4 zz^ l/ tr s/c [5Q(z,0)T37 5 ^Q(0,z)7 5 T 3 ]. 

(IV.7) 

Similar expressions for Q^ v and T^" are also presented 
in 0, [34[ for a single-flavor NJL model. To study 
the effect of very strong magnetic fields, the authors 
of 0, HiJ use the fermion propagator, arising from 
Schwinger proper-time method [351 ] . in the LLL ap- 
proximation. In the present paper, however, we are 
interested on the full eB dependence of these coeffi- 
cients for the whole range of eB £ [0, 1] GeV 2 , and 
have to consider, in contrast to @, Hj], the contribu- 
tions of higher Landau levels too. To do this, we use 
the Ritus fermion propagator 



Sq{x,v) = / V P e 

p=0 



-ip-(x-y) 



xP^D^ip) P p ( yi ), (IV.8) 

arising from the solution of Dirac equation in the 
presence of uniform magnetic field using Ritus eigen- 
function method. The same expression for Sq(x,u) 
appears also in [36| . In (IIV.8|) . p = (p , 0,p 2 ,p 3 ), 
Vp ee , and P p ( Xl ) is given by 

P P (xi) = ll.f+ s (xi) + Il P f- s (xi)] 

+f [/ P +s (^i)-n P /- s (x 1 )] 7 V,(iv.9) 

where, s = sgn(Qe-B), and n p ee 1 — <5 p0 considers the 
spin degeneracy in the LLL. The functions fp S (x\) 
are defined by 

f+'fa) = 4> p (x 1 -sp 2 £%), p = 0,1,2,..-, 
fp s ( x i) = <f) P -i (x! - sp 2 £ 2 B ) , p = 1,2,3,- ••, 

(IV.10) 

where <fr p (x) is a function of Hermite polynomials 
Hp(x) in the form 



4> p (x) = dp exp 



2£% 



(IV.ll) 



Here, a p ee (2 p p\^/tt£b)^ 1 ^ 2 is the normalization fac- 
tor and Ib = \QeB\~ x l 2 is the magnetic length. In 
(|IV.8p . Dq(p) ee 7 • pQ — m, with the Ritus four- 
momentum from (jlll.llj) . Note that since Q is a 2 x 2 
matrix in the flavor space, and therefore P p (xi) 
are matrices in the flavor space. In what follows, we 
will first determine (M%,M% ) and (Q^, T^ v ) at zero 
(T, \i) and in the presence of a constant magnetic field. 



We then introduce T and /i using standard replace- 
ments 



Po = i(2n + 1)ttT — //, and / 



dpp 
2ir 



n 

(IV.12) 



and present the result for (M 2 ,M 2 ) and {Q» v ,T* V ) 
at finite (T, /i, eB) up to an integration over p^- 
momentum and a summation over Landau levels p. 



A. (M%,M* ) at finite (T,fj,,eB) 
1. Ml at finite (T,/x,eB) 



To compute M 2 from pV.4l) . we use the definition of 
the Ritus fermion propagator PV.8I) . and arrive first 
at 



M. 



j2 



l 

2G 



/ d'z J2 V P Vk 



-iz-(p—k) 



•,k=Q 



xtr sc (D- 1 {p)Pp(Q)K k {Q)D-\k)K k {z 1 )Pp{z l )) . 

(IV.13) 

Here, the summation over q £ {§> — g} replaces the 
trace in the flavor space, and in D q , q are the eigen- 
values of the charge matrix Q = diag(2/3, — 1/3). 
After performing the integration over = 0,2,3, 
and using the definition of D^ 1 as well as the Ritus- 
momentum (jlll.lip . with Q replaced by q, we get 



Mi 



2G 
x / dp 2 tx 



q p,k=0 y ' 



1-Pq 
7(0) 



-I P k(p2,k 2 ) 



7 • k q — m 



XJ£'(*2,P2) 



(IV.14) 



k—p 



where the factor 3 behind the integral arises from the 
trace in the color space using tr c (Ij\r eX Ar c ) = 3, and 
two functions I p k and jf® in (|IV.14|) are given by 

I P k(p 2 ,k 2 ) ee P p (0)K k (0), 
J { kp \k 2 ,p 2 ) ee f dziK k { Zl )P p {zx). (IV.15) 



Here, K k (x\) is defined similar to P p (xi) from pV.9() 
Kk{xi) = i[^ s (^i) + n fe .g,T s (xi)] 

+|[^ s (^)-n fc5 r(^)hV,(iv.i6) 

with g^ s (xi) defined as in (|IV.10[) . withp 2 replaced by 
k 2 . Note that for k 2 — p 2 , which is included in the con- 
dition k = p in PV.14|) , we have g k s \k 2 =p 2 = f k S - I* 1 
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what follows, we will first evaluate the integration over 
Z\ in (|IV.15|) . Using then the orthonormality of the 
Hcrmite polynomials appearing in P p from (|IV.9[) , the 
^-integration can also be performed. We will eventu- 
ally end with an expression for M 2 , that includes only 
two integrations over po and p 3 momenta. To start, 
let us first rewrite I p k and jf® from PV.15I) using the 
definition of P p (xi) from pV\9|) and pV.16|) . We get 

I p k(p2,k 2 ) = a+ k (p 2 ,k 2 ) + isj 1 j 2 a~ k (p 2 ,k 2 ), 



4°J(k 2 ,p 2 ) ee Al^(k 2 ,p 2 )+is^A-^(k 2 , P2 ), 

(IV.17) 



-(o) 



-(0)/ 



Here, p 2 q = p\ — 2\qeB\p — p 2 and for p = k, k 2 = 
Pq — 2\qeB\k — p 2 . To perform the integration over p 2 , 
we first compute 



W. 



(0) 



■pk 



^ 2 / p +s (0)/+ s (0). 



(IV.23) 



This can be done using the definition of / p l_s (0) 
in terms of Hcrmite polynomials [see (|IV.10[) and 
PV.IIP ]. and their orthonormality relation 



d£ e~ e H p (e)H k (£) = 



Opk 



(IV.24) 



where 



<*%{P2M = -[/+ s (o).g+ s (o)±H p n fc /- s (o) 5 r(o)] 
Ati°\k 2 , P2 ) 



- / dz 1 [f+ s (z 1 )g+ s (z 1 )±n p Tl k f p s (z 1 )g k - s (z 1 )]. 

(IV.18) 



In this way, the integration over z\ in pV.15[) reduces 
to an integration over z\ in A kp ^° (k 2l p 2 ). The latter 
can be performed using 



dz 1 f+ s (z 1 )g+ s (z 1 ) 



— p„ — a 



(-l)P2 k a k - p 

y/2 k +Pk\p\ 



-U (-p,l + k-p, 2a 2 ) ,(IV.19) 



where a = tBlP2 n ~ K2j and l B = \qeB\~ l l 2 , and 



_ lg(P2-fc2) 

— 2 

U(m, n, z) is the confluent hypergcomctric function of 
the second kind [33] ■ This can, however, be simplified 
by implementing the condition k 2 = p 2 , which is re- 
quired in pV.14l) . In this case a vanishes, and PV.19P 
therefore reduces to 



dZ 1 f+ S (z 1 ) 9 + S (z 1 ) 



k 2 =P2 



= S pk . (IV.20) 



Plugging this result in PV.18P and using H p = n p , we 
arrive at 



P2) 



-(l±U p )6 pk . 



(IV.21) 



Plugging further PV.17I) in pV.14[) . and performing 
the traces over the 7-matrices, using tr s (7 M 7 t/ ) = Ag^ 

and tr s {riiSlvlpla) = 4 (fiVSV ~ 9 PP 9ua + 9^9u P ), 
the er-meson squared mass matrix is given by 



Mi = — + I2i 



2G 



dp 2 



q p,k=0 



(« 0) + 



dp dp 3 
(2tt) 3 



a pk^kp 



(m 2 +p q ■ k q ) 



{Pi 



2 W q 



m 2 ) 



leading to 



TT", 



(0) 



pk 



2fe-l)P+* / dp' 2 e-^HM)H k (p' 2 ) 



Ib 
5 pk 

1 B 



(IV.25) 



with p' 2 = £bP2- Moreover, we arrive at the useful 
relation 



d Pz a % fa M)A^ p 0) {k 2 ,p 2 )\ k 
5 P k 



Af 



a ( 1±n p) 



2=P2 

2 



(IV.26) 



arising from pV.23p . Plugging these results in (IIV.22|) 
and summing over k, the cr-meson squared mass ma- 
trix at zero temperature, chemical potential and non- 
vanishing magnetic field is given by 



9 1 
M = — 
2G 



-6» e \i eB \J2 a p [ 



dp a dp 3 (p 2 q + m 2 



) 3 (Pq ~ TTi 2 ) 2 ' 

(IV.27) 



where a p = 1 + n p is the same spin degeneracy factor 
that appears in pil,14[) . To introduce the temperature 
T and the chemical potential fi, we use the method 
described at the beginning of this section [see pV.12j) ]. 
The mass squared matrix corresponding to er-meson 
at finite (T, fx, eB) is therefore given by 



Ml = 2^-6 E \<ieB\ 



9£{§.-|} 



x f> / 72 %2^ 1 (0) K) + 2m 2 5f (o; p )],(IV.28) 

where w 2 =p 2 + 2\qeB\p + m 2 , and sf^ (u) p ), l=\,2 
are defined by 



2p 2 k 2 a k A 



-(0) 



pk kp 



(P 2 q 



2 )(k 2 



m 2 ) 



(IV.22) 



p—k 



Si ^ (Mr, 



* E 



n— — 00 



Oil 



2\m 



(p 2 o-" 2 p y 



(IV.29) 
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with £ > 1, m > 0. Using 



Sj 0) K) = ^-[l-itywp)], (IV.30) 

and assuming that = 0, Vm > 0, following recur- 
sion relations can be used to evaluate S^ m \ui p ) from 
PV.29|) for all £ > 1 and m > 0, 



1 



sj%p) 



2(^ — l)u> p duip 



W > 2, 



5, (m) K) = 5^r 1) K)+^ (m - 1) K).(IV.31) 

In pV.30|) . Nf(u p ) = n~f(Lu p ) + nj(io p ) and n^(uj p ) 
are fermionic distribution functions 



(IV.32) 



In the following paragraph, the same method will be 
used to determine M%„ at zero and nonzero (T, /i) and 
for non- vanishing eB. 



2. Mlo at finite (T, fi, eB) 

To determine the squared mass matrix Mj from 
(|IV.5p . corresponding to 7r°, we use the definition of 
the fermion propagator PV.8I) - PV.9|) . and arrive first 
at 



= ^3 - * / d * z E I V P Vk e~ iz ' [§ ~ h) 
D Q 1 (p)Pp(0)T 3l5 K k (0)D Q 1 Ck) 



2G 

X tTgfc 



xK k (zi)j 5 T 3 P p (zi) 



(IV.33) 



Using the anticommutation relation {75, 7 M } = lead- 
ing to [75,-Kfe] = 0, we simplify first the combination 
j5K k (0)Dp 1 (k)K k (zi)j 5 in (IIV.33|) . and arrive at 



75^ fc (0)i? Q 1 (fc)^(z 1 )75 = -K k (0) 



1 



7 • k + m 



K k (zi). 
(IV.34) 



Plugging this relation in PV.33P and performing the 
integration over Zi, i = 0, 2, 3, we arrive at 



^o = - + 3, 



E 

p,k=0 



dpadp 3 . 

J- I d P2 tv fs 



(2tt) 2 



7 • p — m 



-I p k(p2,k 2 ) 



XT 3 



1 



7 • k + m 



T 3 J k p\k 2 ,P2) 



(IV.35) 



k—p 



where I pk (p2, k%) and (£2,^2) are given in pV.15|) . 
We follow the same method leading from pV.14|) to 



PV.27P to evaluate the traces over the 7-matrices and 
to perform the integrations over z\ and p 2 in PV.35P . 
We arrive after a lengthy but straightforward compu- 
tation at 



2G 



6i ^ \qeB\ 



OO 



dpodps 



(IV.36) 



Thus, the mass squared matrix corresponding to 7r° 
at finite (T, /j,, e-B) is given by 



M *° = E 1^1 

3£{| -|} 



x E«p/7^^ 0) K)' < TV -'-> 7 » 

p=o J [Z7T> 



where s[°\u; p ) is given in pV.30p . In Sec. [Vj the 
integration over p% and the summation over Landau 
level p, appearing in pV.37[) . will be performed nu- 
merically. 



B. (g^,T^) at finite (T,fi,eB) 
1. at finite (T,fi,eB) 



We start by computing from PV.6P at zero (T, /z) 
but non-vanishing eB. To do this, we use the defini- 
tion of the Ritus propagator pV.8p , and arrive first 
at 



Q>lV = 9 E / d ^ z ^ E / V P ™ e -*«-(P-*) 



p,k=0 



xtr s 



D- 1 (p)P p (0)K k (0)D- 1 (k)K k (z 1 )P p (z 1 ) 



(IV.38) 



After performing the integration over Zi,i = 0,2,3, 
and using the definition of D~ x , the diagonal elements 
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of Q^ v are given by 



q k,r=0 ' 



8Q 



(7 • r q - to) 



1 



EE 

q k,p=0 



Xlpk(P2,k 2 ) 



(7 • kg - TO,) 

2?p tr 

1 



(l-Pq 



dpi 



EE 

q k,p=0 
Ipk(j>2,k 2 ) 



(7 • k q - TO) 

T>p tr 



4p(k 2 , P2 ) 



k—p 



1 



(7 • fc 8 - m) 



(7 • p q - to) 
4°\k 2 , P2 ) 



k—p 

(IV.39) 



In j — 0, 3, r g = k q +£ q and r = k + £. Moreover, 
two functions I p k(p 2 ,k 2 ) and jj^ p \k 2 ,p 2 ) are defined 
in (jlV.lSp . and J^ p \k 2 ,p 2 ) is denned by 

4 2 J(k 2)P2 )= I ' dzxzlK^Ppizx), (IV.40) 



with P p and fTfc given in (IIV.9|) and PV.16|) . respec- 
tively. Following the method presented in the first 
part of this section, leading from (|IV.14|) to (|IV.27[) , all 
non-diagonal elements of Q^ u turn out to vanish, and 
therefore, as it is claimed in Sec. HH Q^ v = Q^g^ v (no 
summation over fj,). This is similar to what also hap- 
pens in the single-flavor NJL model [2j. We therefore 
focus on gw,n = 0, • • • , 3 from (IIV.39|) . which shall 
be evaluated using the same method as before. Eval- 
uating the fc-integration in , j = 0, 3 from (jIV.39|) . 
using an additional Feynman parametrization, we ar- 
rive first at 

gOO = _ g 33 = 3. £ \qeB\ 

9e{f,-|} 
dp dp 3 f 1 



E< 

p=0 



(27r) 3 |_ — TO 2 ) 2 3 (pi — TO 2 ) 3 



(IV.41) 



At finite (T, /i, eB), we therefore have 

s 00 = -a 33 = -3 J2 \* eB \ 

96{|.-i} 

(IV.42) 

To determine G 11 from (|IV.39j) . we shall first evaluate 
the Z\ integration in jj^ p \k 2 ,p 2 ) from (|IV.40|) at k 2 = 



p 2l as it is required from pV.39|) . To do this, we first 
define 

4 2 J(k 2 , P2 ) ee At p {2 \k 2 , P2 ) + ll 1 1 ~sA- p {2 \k 2 , P2 ), 



with 



A^ 2 \k 2 , P2 ) 



(IV.43) 



(IV.44) 



= o £k P {k 2 ,p 2 ) ± U p Ilk£k~i,p-i(k 2 ,p 2 ) 

Here, L k p(k 2 ,p 2 ) is defined by 

£ kp (k 2 , P2 ) = J dzt [z l5 + s (zi)][zi/+ s (zi)].(IV.45) 

To determine C kp for k 2 = p 2 , we use the definition of 
fp S from (IIV.10I) in terms of the Hermite polynomials 

Hp and their standard recursion relations dH ^ x ^ = 
2kHk_i(x) and H k +i(x) — 2xH k (x) — 2kHk-i(x), to 
arrive first at 

= t B (Cr+if+tM) + c P fp-M) +v' 2 f; s {zi)) , 

(IV.46) 

where C p = \ and p 2 = £bP2- Replacing (|IV.46[) in 
(|IV.45[) . setting k 2 = p 2 , and integrating over z%, we 
get 



£k P (k 2 = p 2 ,p 2 ) = l\ 



[ptp+l +p7) hp 



+CpC p -\5k,p~2 + C p +iCp+ 2 Sk t p+ 2 + 2p' 2 (C p 5k lP - 



+Cp+i5k,p+i) 



Thus, A^ 2 \p 2 ,p 2 ) in (jIV.431) are given by 



(IV.47) 



A±(2 



hp 

(k 2 =p 2 ,p 2 ) = -y 



C ± S. 



kp 



+C p _i (Cp ± lJpUkC p - 2 ) 5k, P -2 
+C p +i (Cp+2 ± flpilfeCp) Sk, p + 2 
+2p' 2 (C p ± npiT fc Cp_i) 5k, p -i 



-2p 2 (Cp+i ± ripilfcCp) 6k,p+i 



,(IV.48) 



where the coefficients C^ = +p 2 (l ± Hp) with 
D ± ee Cl p+1 ± Cf^IIp. Plugging pV^8| in g glj 
and the resulting expression in Q 11 from (|IV.39|i . and 
performing the trace over 7-matrices, we arrive at 

s = 6 *EE J 

q p,k=0 J y ' 



dp2 



(< k A + k i 2) +a; k A-^)( P q.k q+m ^ 



(Pq — TO 2 )(fc 2 — TO 2 ) 



+ 



2p 2 k 2 a pk A kp 



(2) 



(p 2 — TO 2 )(fc 2 — TO 2 ) 



(IV.49) 



k—p 



11 



where a^ k are defined in PV.18P . The integration over 
P2 is then performed using 



(i) _ 



pk 



(2) 



d P 2P' 2 f+ s (0)f+ s (0) 
1 



i 2 

B 



W/ P +s (o)/ fc +s (o) 



— +72" (Clp+l^P + C p+ 2C p+ i8k.p+2 



+C p Cp-i5k,p-2) 



(IV.50) 



These results arise from the orthonormality relations 
of the Hermite polynomials (|IV.24[) , in the same way 
that Wj£> from pV.23|) is derived. Using W$ and 



wf k ] from PV.50I) . we get 

J dp 2 a± k A k £ 2) \ k= . = ^ 



D ± + C| , j ± ITpCjp.! 



x (i ± n p )4 P - 2 (c p ± UpUkCp^y 4,p-i 



-2 (Cp + i ± UpILkCp) 2 5k,p+\ 



(IV.51) 



Plugging these relations in (|IV.49[) . we finally arrive 
at 



Zi 

f dp Q dp 3 f (p q ■ k q + m 2 ) (1) 

X Z^2^ / (2tt) 3 \(p2_ m 2)(fc2_ m 2) Pfc 



9 p,fc=0' 

p 2 k 2 



(p 2 — m 2 )(fc 2 — m 



(2) 



2\ w pfc 



(IV.52) 



k—p 



where 

C# = [(2p+l)+IL p (2p-l)]S kp 

-[p + n p n fe (p- i)]4 iP _i - [(p+ 1) + u p u. k p]dk,p+i, 



C { $ = -[(2p + 1) - (2p - l)n p ](l - 



-b + n p n fc ( P - i - 2Vp(p - i))]*fc, P -i 



+\{p + i) + n p n fe (p - 2 a /p(p + i))]4, p +i. (iv.53) 

To determine £/ 22 from pV.39[) . we perform the traces 
over the 7-matrices and arrive first at 

dp dp 3 



->22 



(2tt) 3 



x / dp 2 



6i £ £ 

96{ |,_i } p,fc=0' 

<'^ ? .fc g + m 2 )ivW+2p 2 fc 2 ^ ) 



(p^ — m 2 )(k 2 — m 2 ) 



(IV.54) 



where 



d 2 , 

d^A a 



a: 



-(o) 



dp 2 

d 2 

^1 l«p fc ^p 



^pk^kp 



&kA- k r) 



Plugging the definitions of ct^ k (p 2 ,k 2 ) and 



A k ^ 0) (p2,k 2 ) from (IrVTSll in (jlVlIji . and per- 



forming the integration over p 2 in pV.54[) by making 
use of from pV.251) . we arrive after a lengthy 

but straightforward computation at 

J dp 2 N$(p 2 ,k 2 =p 2 ) = ~\cf k \ 



dp 2 N^ k \ P2 ,k 2 =P 2 ) = iC'l', (IV.55) 



(2) 



pk 



4 P k 



where C& and C p 2 k are given in pV.53j) . This leads 
eventually to 

(IV.56) 



G 22 =9 1 \ 



with Q given in PV.52|) . Note that the equality 
Q 11 = Q 22 arises also in a single-flavor NJL model 
in 0, where the form factors of the effective kinetic 
term are computed at zero temperature and chemi- 
cal potential and in the regime of LLL dominance. 
In pV.56j) . this regime is characterized by k = p = 0, 
where k and p label the Landau levels. At finite (T, 
Q 11 = Q 22 is therefore given by 

00 „ , 

g" = g 22 = -3 £ E/t^{^ 0) K) 



13 ' 3 J 

1)2 _2 , „1 C (0), 



+2[(2p+ l)4m 2 +p] S^K) 
-2[(2p+l)4m 2 + (p+l)] ^% p+1 ) 
+<S p0 [<S} 0) + 2m 2 ^ 0) (u p )] } , (IV.57) 

where S{° (cj p ) is given in pV.30|) and <sij (wp) can 
be evaluated using the recursion relations (IIV.31|) . 

2. J rM " at finite (T, n, eB) 

We start the computation of the elements of the ma- 
trix T^" by considering its definition from pV.7[) . and 
arrive after plugging the Ritus propagator pV.8[) in 

dTVTrt at 



I 

xtr 



d z z 



p,k=0 



sfc 



DQ 1 (p)P p (0)T 3l5 K k (0)DQ 1 (k) 



x K k (zi)j 5 T 3 P p (zi) 



(IV.58) 



Using pV.34[) and following the same method as is 
used to determine Q^ v in the previous section, we ar- 
rive after some work at 

j-oo = _^33 = 3l \qeB\ 

{f.-i> 

dp dp 3 1 



OO 



£• 

p=0 



(2tt) 3 (p 2 -m 2 ) 2 ' 



(IV.59) 
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and 



T 11 = F 22 



3*EE 

q p,k=o 



dp dp3 



(2tt) 3 
2\qeB\y/pJe 



iPg ■ k q - m 2 ) (1) 



(p2- m 2)(fc2_ m 2) Pk 



(p 2 - m 2 )(k 2 - m 2 ) pfe 



(2) 



(IV.60) 



k—p 



At finite (T, J 7 ^, ^ = 0, • • • ,3 are therefore given 

by 



r 



■00 



-T 



:S:S 



-3 J] keB\±a p [^S^^ p ), 



(IV.61) 



as well as 

J" 11 = J" 22 

OO 

—3 E E 

:(o) 



d P* J or.. o(o). 



(27T) 



2[p5| U) (u; F 



-(p + 1) Sl u > p+ i)] +<5 P o5i U) (^)|.(IV.62) 

All non-diagonal elements of T^ v turn out to van- 
ish. As we have described before, the remaining p$- 
integration and the summation over Landau levels ap- 
pearing in the final results of Sees. II V Al and II V Bl for 
the squared mass matrices (M 2 , M 2 ) as well as form 
factors (kinetic coefficients) {Q^ v , F^), will be eval- 
uated numerically in the next section. Using these 
results, the (T, fi, eB) dependence of pole and screen- 
ing masses as well as the refraction indices of neutral 
mesons will be explored. 



V. NUMERICAL RESULTS 

In Sec. IIIII we have introduced the one-loop effective 
action r off [a, 7?] of a two- flavor NJL model describing 
the dynamics of non-interacting a and tt mesons in 
a hot and magnetized medium. We have then deter- 
mined the corresponding one-loop effective potential 
of this model il eff (m; T, /z, eB), up to an integration 
over p3-momentum and a summation over Landau lev- 
els, labeled by p. According to our description in Sec. 
HI1 the global minima of Q e , f (m; T, /i, eB) can be used 
to determine the squared mass matrices (M 2 ,M 2 ) 
and the coefficients of the form factors (kinetic coef- 
ficients) (G^, F^ v ), corresponding to neutral mesons 
and appearing in the effective action (lIV.ip . In Sec. 
IIVI we have described the analytical method leading 
to (M 2 ,M 2 ) and [G^ ,F» V ) at finit e {T^ eB). The 
squared mass matrices are given in (|IV.28[) as well as 
(|IV 37j) an d the form factors in (|IV.57|i . (|IV.6ip as well 
as (|IV.62|) . All these results are presented up to an 



integration over p3-momentum and a summation over 
Landau levels p. In this section, we will first use the 
one-loop effective potential f|III. L4[) . to determine nu- 
merically the (T, /x, ei?)-dependence of the constituent 
quark mass m = mo + co for non- vanishing bare quark 
mass mo- This will be done in Sec. IV Al by keeping 
one of these three parameters fixed and varying two 
other parameters. We then continue to explore the 
complete phase portrait of our magnetized two-flavor 
NJL model in the chiral limit mo — > 0. Our results are 
comparable with the results previously presented in 
[39l lioj . Similar results are also obtained in || , where 
the two-flavor NJL model, used in the present paper, 
is considered with additional diquark degrees of free- 
dom to study the chiral and color-superconductivity 
phases in a hot and magnetized quark matter. In 
Sec. IV Bl we will then evaluate the above mentioned 
P3-integration and the summation over Landau levels 
numerically. This gives us the possibility to study, in 
particular, the T-dependence of (M 2 , M 2 ) as well as 
(g^, T» v ) for fi = and various eB = 0, 0.03, 0.2, 0.3 
GeV 2 (or equivalently eB ~ 0, 1.5m 2 , 10.5m 2 , 15.7m 2 
for m,,- = 138 MeV). As we have described in Sec. HI 
the magnetic fields produced in the con-central heavy 
ion collisions at RHIC and LHC are estimated to be 
in the order of eB ~ 1.5m 2 and eB ~ 15m 2 (or equiv- 
alently, eB ~ 0.03 GeV 2 and eB - 0.3 GeV 2 , respec- 
tively) num. 

Hence, our results for small values of 
magnetic fields (here, eB = 0.03 GeV 2 ) may be rele- 
vant for the physics of heavy ion collisions at RHIC, 
while our results in the intermediate magnetic fields 
(here, eB = 0.2, 0.3 GeV 2 ) seem to be relevant for the 
heavy ion collision at LHC. In Sec. IV CI we will finally 
present a number of applications of the results pre- 
sented in the second part of this section. In particular, 
we will determine the T-dependence of the pole mass 
as well as the refraction index and screening mass of 
neutral mesons for fi = and eB = 0,0.03,0.2,0.3 
GeV 2 . To do this, we will use the corresponding dis- 
persion relations of a- and 7r° mesons. The goal is to 
study the effect of uniform magnetic fields on meson 
masses and refraction indices and explore the inter- 
play between the effects of temperature and the ex- 
ternal magnetic fields on these quantities. We will, in 
particular, show that uniform magnetic fields induce a 
certain anisotropy in the mesons refraction indices and 
the screening masses in the longitudinal and trans- 
verse directions with respect to the external magnetic 
field. Detailed studies on eB and fi dependence of all 
the above physical quantities, together with other pos- 
sible applications of (M 2 ,M 2 ) and {Q» v , F^), e.g. 
in studying the mass splitting between charged pion 
masses will be presented elsewhere (30j. 7 



7 The mass splitting between 7r+ and tt~ is recently discussed 
in [20I , l2l| , using chiral perturbation theory in the presence 
of constant magnetic field. 
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A. Chiral condensate and complete phase 
portrait of a magnetized and hot two-flavor NJL 
model in the chiral limit 



Using the thermodynamic potential from (jlll. 14[) . we 
will determine, in what follows, the chiral condensate 
and the complete phase portrait of the two-flavor NJL 
model at finite T, \x and eB. The notations and math- 
ematical method used in this paragraph are similar 
to what was previously used in Q. To determine the 
chiral condensate, we have to solve the gap equation 
numerically 



dn ett (fh; T, fz, eB) 



dm 



= 0. 



(V.l) 



Here, m — rriQ + a and m — mo + 0o, as are introduced 
in Sec. IIII1 Our specific choice of parameters is [4l[ 

A = 0.6643 GcV, G = 4.668 GeV _2 ,m = 5 MeV, 

(V.2) 

where A is the UV momentum cutoff and G is the NJL 
(chiral) coupling constant. To perform the momentum 
integration over p and j»3, we have introduced, as in 
81, smooth cutoff functions 



1 



fP 



l + ex P (M=A) 
1 



>A,B 



1 + cxp( VpI+M^-A 



(V.3) 



corresponding to integrals with vanishing and non- 
vanishing magnetic fields, respectively. In p 
labels the Landau levels. Moreover, A is a free pa- 
rameter, which determines the sharpness of the cutoff 
scheme. It is chosen to be A = 0.05A, with A given 
in (|V.2[) . Using the above smooth cutoff procedure, 
the above choice of parameters leads for vanishing 
magnetic field and at T = fj, = to the constituent 
mass m ~ 308 MeV. 8 Let us notice that the solutions 
of (|V.1[) are in general "local" minima of the theory. 
Keeping <tq ^ and looking for "global" minima of the 
system described by fi off (m; T, fi, eB) from piI.14[) . it 
turns out that only in the regime /i G [0, 350] MeV, 
T G [0,390] MeV and eB G [0,0.8] GeV 2 , the global 
minima of f2 ef{ are described by nonzero Co- In these 
regimes, the chiral symmetry is spontaneously broken 
by non- vanishing ctq. 9 All our numerical computa- 
tions in the present section are therefore limited to 
these regimes. Note that because of non-vanishing 



quark mass mo, the transition from the chiral sym- 
metry broken phase to the normal phase is a smooth 
crossover (see the descriptions below). 

In Fig. [TJ the T, /i and eB dependence of m are 
presented. In Fig. HJa), the T-dependence of m is 
demonstrated for fixed /i = and eB = 0,0.2,0.5 
GeV 2 . Although the transition from the chiral sym- 
metry broken phase, with m ^ to the normal phase, 
with m ~ mo ~ 0, is a smooth crossover, but as 
it turns out, for stronger magnetic fields the tran- 
sition to the normal phase occurs for larger values 
of T, whereas for eB = 0, this transition tempera- 
ture into the crossover region is smaller. Moreover, 
at T G [0, 100] MeV, where m is almost constant, the 
value of m increases with increasing eB. All these 
effects are related with the phenomena of magnetic 
catalysis 0, @], according to which, magnetic fields 
enhance the production of cto ~ (ipip) condensate, 
even for very small coupling between the fermions, 
and therefore catalyze the dynamical chiral symme- 
try breaking. Similar effects occur also in Fig. [TJb), 
where m is plotted as a function of fi, at fixed T = 120 
MeV and for various eB = 0, 0.2, 0.5 GeV 2 . At fi = 0, 
for instance, the value of m increases with increasing 
eB. In Fig. HJc), the e-B-dependence of m is demon- 
strated for fixed \i = and T = 60, 180 and 220 MeV. 
As it turns out, for fixed value of eB, m decreases 
with increasing T, and as it turns out, this "melt- 
ing" effect persists in the whole range of eB G [0,0.8] 
GeV 2 , although it is partly compensated by the mag- 
netic field in the regime eB > 0.6 GeV 2 . Let us notice 
that, according to our results in 0, Q, for a certain 
threshold magnetic field eB t ~ 0.45 GeV 2 , the mag- 
netic field is strong enough and forces the dynamics 
of the system to be mainly described by the LLL. In 
this regime, m increases linearly with increasing eB 
[see Fig. HJc)]. Later, in [27| . the threshold magnetic 
field is estimated to be in the order of B ~ 10 19 Gaufi. 
In the present paper, however, the threshold magnetic 
field turns out to be eB t > 0.7 GeV 2 [or equivalently 
B ~ 1.2 X 10 20 Gaufi]. 10 The T and eB dependence of 
m at fixed chemical potential \i and various eB and T 
are discussed recently in [38| using lattice gauge theory 
methods in the presence of constant (electro)magnetic 
fields. Our original results from Q as well as the re- 
sults presented in Figs. [UJa) and [He) are consistent 
with the results arising from lattice simulations [38j . 

According to our results in Q, in the chiral limit 
mo — > and for vanishing magnetic field, at high tem- 
perature and small chemical potential, the transition 
from the chiral symmetry broken to the normal phase 
is of second order. In contrast, at low temperatures 



8 For sharp UV-cutoff, m turns out to be m ~ 300 MeV, as 
expected. 

9 For mo 7^ the chiral symmetry of the original Lagrangian 
is explicitly broken. 



The exact value of threshold magnetic field eBt is determined 
from [ r^g, J = 0, where [a] is the greatest integer less than 

or equal to a. For up quark eBt — 0.67 GeV 2 and for down 
quark eB t ~ 1.33 GeV 2 . 
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FIG. 1: (a) and (b) The T and /i dependence of the constituent quark mass m = mo + uo for fixed eB — 0, 0.2, 0.5 GeV 2 
and for fixed \x — and T = 120, respectively. Here, mo — 5 MeV is the bare quark mass and cto is the chiral condensate, 
(c) The eB dependence of m is demonstrated for fixed /i = and various T = 60, 180 and 220 MeV. 
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FIG. 2: Complete phase portrait of a two-flavor magnetized NJL model at finite T, [i and eB in the chiral limit of 
vanishing quark mass mo. Blue solid (green dashed) lines denote the second (first) phase order phase transition between 
the chiral symmetry broken and normal phases. Two branches of the first order critical line for fj, — 320 MeV are denoted 
by green double-dashed lines. 



and higher densities, the second order phase transi- 
tion goes over to a first order one. In the presence 
of a uniform magnetic field, this picture remains es- 
sentially the same. The only difference is that for 
/i = 0, the transition temperature increases with in- 
creasing eB. Moreover, for eB ^ 0, the second or- 
der phase transition occurs at higher temperatures 
and lower densities comparing to the case of vanish- 
ing magnetic fields. These two effects of the uniform 
magnetic field on the T — \x phase diagram of a two- 
flavor NJL model in the chiral limit are demonstrated 
in Fig. HJa). Both effects are manifestations of the 
phenomenon of magnetic catalysis in the presence of 
constant magnetic fields [l|, |2|. In all the plots of Fig. 
[2 the green dashed (blue solid) lines denote first (sec- 
ond) order phase transitions. To determine the first 
and second order phase transitions, the method de- 
scribed in 0, d, [39l l4fj( is used. The first order critical 
lines between the chiral symmetry breaking and the 
normal phase is determined by solving 



9f2 eff (m; T, n, eB) 



dm 



0, 



simultaneously 11 The second order critical line be- 
tween these two phases is determined using 



lim 



dri Btt (m; T, n, eB) 
dm 2 



0. 



(V.5) 



To make sure that after the second order phase tran- 
sition the global minima of the effective potential are 
shifted to m = in (|V.5[) . and in order to avoid insta- 
bilities, an analysis similar to [42| is also performed. 

In Fig. HJb), the T — eB phase diagram of our 
model is plotted for various fi = 240,280,320,340 
MeV. Let us notice that for (j, = 320 MeV [dashed- 
dotted lines in Fig. Efb)], the first order critical 
line has two branches - the first one for eB < 0.1 
GeV 2 and the second one for eB > 0.5 GeV 2 , at 
relatively low temperature. In the intermediate re- 
gion 0.1 < eB < 0.5 GeV 2 , the chiral symmetry 
breaking phase is disfavored. In Q, we have stud- 
ied the T — eB phase diagram of a two-flavor NJL 
model including meson and diquark condensates. We 
have shown that in the above mentioned intermediate 
regime 0.1 < eB < 0.5 GeV 2 at low temperature and 
for /i = 320 MeV, the two-flavor color superconduct- 
ing (2SC) phase is favored. For fj, > 320 MeV, the 



a ff (m + 0; T, fx, eB) = Q eff (m = 0; T, /*, eB), (V.4) 



In the chiral limit mo — > 0, m \ 
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first branch appearing for /i = 320 MeV and eB < 0.1 
GeV 2 disappears, and the whole region of eB < 0.6 
GeV 2 is favored by either the normal phase, when no 
diquarks exist in the model, or by the 2SC supercon- 
ducting phase, when the model includes both meson 
and diquark condensates (see Fig. 14 of [8]). 

In Fig. [UJc) , the /i — eB phase diagram of our two- 
flavor NJL model including chiral condensates (a, 7?) 
is plotted for various T = 30, 60, 100, 150 MeV. At 
very low temperature, T < 100 MeV, the transition 
between the chiral symmetry breaking and normal 
phase is of first order (green dashed lines). Whereas 
at these temperatures and for eB < 0.1 GeV 2 , the 
critical /i is almost constant, it decreases by increas- 
ing the strength of the magnetic field in the regime 

0. 1 < eB < 0.4 GeV 2 . This effect, which is for the 
first time observed in 0, [|[ , and later also in [27j , is 
called the "inverse magnetic catalysis", according to 
which at low temperature, the addition of the mag- 
netic field decreases the critical chemical potential for 
chiral symmetry restoration [1, H3|- However, by in- 
creasing the magnetic field up to eB > 0.5 GeV 2 , 

1. e. by entering the regime of LLL dominance, this 
effect is disfavored, so that /i c again increases with 
increasing the strength of the magnetic field. Let us 
also note that similar phenomenon of inverse magnetic 
catalysis appears also in Fig. [2jb), where for fixed 
/i = 280 MeV, the first order critical line T c (green 
dashed line between C\ and C2) decreases with in- 
creasing eB from 0.1 < eB < 0.3 GeV 2 and contin- 
ues to grow up with increasing the strength of the 
magnetic field up to regime of LLL dominance, i.e 
eB > 0.5 GeV 2 . The inverse magnetic catalysis effect 
may be related to the well-known van-alphen-de Haas 
oscillations, which occur whenever Landau levels pass 
the quark Fermi level [43[ . Similar effects are also ob- 
served in @, H, E3|. At higher temperature T > 100 
MeV and for eB smaller than a certain critical eB c , 
there is a second order phase transition between the 
chiral symmetry broken and normal phases (see the 
blue solid lines in Fig. C^c) for T = 100, 150 MeV, that 
replace the green dashed lines for T < 100 MeV). The 
critical magnetic field eB c , for which the second or- 
der phase ends and goes over into a first order phase 
transition is larger for higher temperature [compare 
eB c for two critical points (black bullets) C% and C2 
in Fig. dfc)]. This demonstrates the destructive effect 
of the temperature, which is partly compensated in 
the regime of strong magnetic fields, eB > 0.7 GeV 2 . 
More details on the interplay between three parame- 
ters T, fi and eB on the formation of chiral conden- 
sates ctq in the chiral limit mo —¥ are discussed in 



B. 



(MlM^o) and (Q» v ,T» V ) at 
(T/0,M = 0,eB#0) 



tion, the results presented in (|IV.28[) and (IIV.36I) for 
M 2 and M 2 , and in (1PV421) and (ITV571) for Q» v as 
well as in (IIV.61[) and (|IV62j) for are given up 
to an integration over ^-momentum and a summa- 
tion over Landau levels p. We have performed the p%- 
integration for the set of parameters (A, G, mo) from 
(|V.2[) and the smooth cutoff function (IV. 31) numer- 
ically, and will present the results in what follows. 
In particular, we will present the T-dependence of 
(M 2 ,M 2 ) and (G» v ' ,T» V ) for fixed fi = and var- 
ious eB = 0,0.03,0.2,0.3 GeV 2 . 

Let us start by giving the numerical values of 
(Ml, Ml) and (0^,7^) at T = p. = eB = 0. Ac- 
cording to their definitions in pV.4[l - pV.7[) . where, for 
eB = 0, Sq is to be replaced by the ordinary fermion 

propagator S(z,0) = / (Cj4^^ at zero (T,(i), we 
have the following identities and numerical values 

Ml = 3.656 x 10~ 2 GeV 2 , 
Ml i = 1.734 x 10~ 3 GeV 2 , W = 1, 2, 3, (V.6) 



as well as 



-.01) 



Vi = 1,2,3, 



(T uo ) u = -(r i ) u Vi= 1,2,3, (V.7) 



where 



g QQ = 5.381 x 10~ 2 , 
(J^V = 9.143 x 10~ 2 , W = l,2,3. (V.8) 

Moreover, we have (F ll )n — (T ll ) 2 2 = (J 7lt )33 for all 
£ = (),•■• ,3. 

At finite temperature and vanishing fi and eB, al- 
though the above relations (|V.6I) and (IV7[) between 
different components of Mi as well as Q^ v and T^ v 
are still valid, i.e. we have 



(Mj) u = (M|) 22 = (Af|) 33 , 



(V.9) 



as well as 



As we have described in the first part of this sec- 



goo _ _ —Q 22 = _ g 33 

(T 00 )u = -(T vl )u, V£,i = 1,2,3, (V.10) 

but their values become temperature dependent. In 
Fig. H the T-dependence of (Ml,(Ml) u ) as well 
as (G 00 ,g u ), and ((-F 00 )«, (T u )u) with l,i = 1,2,3 
are plotted for vanishing eB and fi. As it is demon- 
strated in Fig. [3Ja), M 2 and (Mi)u are degenerate at 
T > 220 MeV. This is because the difference between 
these two functions are in terms proportional to the 
constituent quark mass m = mo + 00, that, according 
to Fig. [TJa) almost vanishes in the crossover region 
T > 220 MeV. Later, we will show that the degen- 
eracy of Ml and (M|)« for T > 220 MeV leads to 
the expected degeneracy of a and ir° meson masses 
(m a , m„o) for vanishing eB and fi in the crossover re- 
gion T > 220 MeV 0]. 
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FIG. 3: The T-dependence of (M^,(M|)«) as well as (G 00 ,g"), and ((T 00 )u, (^lee) with i,i = 1,2,3 for vanishing 
magnetic field eB and at fi = 0. 



a) /j = MeV 




o 

_o 
> 

CD 

So 
s 

o 





200 
T (MeV) 

b) pl = MeV 



400 



eB 


= 0.03 GeV 2 


eB 


= 0.2 GeV 2 


eB 


= 0.3 GeV 2 




200 
T (MeV) 



400 



FIG. 4: The coefficient M 2 (panel a) and M 2 (panel b), 
are plotted as functions of T £ [0, 400] MeV at fi — and 
for eB = 0.03, 0.2, 0.3 GeV 2 (red solid, blue dashed and 
black dotted lines, respectively) . 



Let us finally consider the case of (T, eB ^ 0,/z = 0). 12 
As it turns out, the degeneracy in (Mi)u as well as 
g H and (T")«, with i,i = 1, 2, 3 at (T, eB ^ 0, /x = 0) 
breaks down by finite magnetic fields. In other words, 
for (T, eB ^ 0, /i = 0), in contrast to (|V.9[) . we have 



Moreover, in contrast to (|V. 10|) 

£00 = _g33 gll = g22 



(V.12) 



Similarly, in contrast to (|V.10[) . although (-F MM )ii = 



(J-w) 22 ^ (7-^)33, for ah> = 0,- 



(^)« ^ (J" 



, 3, but 

= (T 22 )«, (V.13) 



V i 



and for non- vanishing 
The exact (T, ^, eB) de- 



(M|) n = (M|) 22 ^ (M|) 33 . 



(v.n) 



1,2,3. In Fig. H the T-dependence of M 2 
(panel a) and Af 2 [or equivalently, (Af|) 33 ] (panel 
b) are demonstrated at fi 
eB = 0.03,0.2,0.3 GeV 2 . 

pendence of (M|)h = (M|) 2 2 will be used in [30|, to 
determine the (T, /i, e 5) -dependence of charged pion 
masses. 

In Fig. the T-dependence of G 00 and g 33 
(g 00 = -g 33 ) (panel a) as well as g 11 = g 22 (panel 
b) are plotted for fi = and eB = 0.03, 0.2, 0.3 GeV 2 . 
Whereas g 00 is positive, g 11 = g 22 and C/ 33 are neg- 
ative. Later, we will use the matrix elements of M 2 
from Fig. [Hand the coefficients £7 MM , /i = 0, • • • ,3 from 
Fig. [5j to determine the T-dependence of to ct at [i = 
and for various eB =/= 0. 

In Fig. [5J the T-dependence of , fx = 0, • • • ,3 
[or equivalently, (T^)^} matrices are plotted for 
vanishing chemical potential and eB = 0.03, 0.2, 0.3 
GeV 2 . In the subsequent section, we will in par- 
ticular use M 2 , T 00 and J 733 to determine the T- 
dependence of n° pole and screening masses as well 
as the direction-dependent refraction indices of neu- 
tral pion in the longitudinal and transverse directions 
with respect to the direction of the external magnetic 
field. 



Masses and directional refraction indices of 
neutral mesons 



12 In this paper, we are interested on the effects of magnetic 
fields on the meson masses and their refraction indices at 
T^O and p = 0. The results for T ^ and ^ as well 
as the e-B-dependence of these quantities will be presented 
elsewhere 1 3011 - 



In this section, we will use the results obtained in 
Sec. IV Bl to determine the T-dependence of pole and 
screening masses as well as the direction-dependent 
refraction indices of neutral mesons, a and 7r°, in a 
hot and dense magnetized quark matter. In what fol- 
lows, we will first define these quantities according to 
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FIG. 5: The coefficients Q" 



-a 33 (panel a), Q 1 



(panel b) are plotted as functions of T G [0, 400] MeV for 
H = and eB = 0.03, 0.2, 0.3 GeV 2 . 



the descriptions presented in Sec. |TT] and the corre- 
sponding energy dispersion relations for neutral and 
charged mesons a and 7r^,<? = 1,2,3 mesons [see also 
(flLl~5j) and (IlLTHll ]. 



1 

goo 



a) i 



(^ 00 )i 



■22 



33 



(V.14) 



The pole and screening masses of er-mesons, m a and 



,00 



1,2,3, are then defined by 



Re Ml 



Re G 00 



where, u& , is the directional refraction index of er- 
mesons in the i-th direction, 

1/2 



1/2 



as well as 



,00 - 



„.00 



, (V.15) 



u 



(0 



Re 



Re g 



on 



1,2,3. 



(V.16) 



The (T, /i, ei?)-dependence of m a -,ma and are 
given by plugging M% and g^,^ = 0, • • • , 3 from 
PV.28|) , (|IV.42[) and (|IV.57|) in the above relations. 
As concerns the pions, we choose the basis (7t ± ,tt°) 



instead of the real basis 7? = (tti, 7T2, ^3). Here, 
71"^ = (7Ti ± iTT2)/V2 and 7r° = 7T3. Using this new 
imaginary basis, the energy dispersion relations E 7Tl 



from (fVU4|) for (n- 



turn out to be 



El, 



Elo 



eB(2l+l) 



~ [(^ 00 )iiT*(-F° )i2 

[(^33 )iiTl(J -33 )i2] 



[(M*) llT i(M*U 
[(7-° ) U Ti(^ 00 )i2] 



33 2 
Pi 



" (^" 00 )33 

(^33 2 
(J-00) ^2 



33 



(-^ 33 )3 3 



J.33 



?>3 



(Mg) 33 
(^°°) 



33 



(V.17) 



Note that since ^ are charged pseudoscalar parti- 
cles, their energy dispersion relations in the presence 
of constant magnetic fields have discrete contribu- 
tions. According to our results in (45[, the energy 
levels are labeled by £, in the form given in the first 
term in (|V.17j) . The above dispersion relations for 
charged pions are comparable with the dispersion re- 
lations presented recently in [2l| (see Eq. (2.10) in 
\2l\). According to the formalism presented originally 
in [2J and generalized to a multi-flavor system in the 
present paper, in contrast to the relations presented 
in [2l[ for charged pions, the nontrivial form factors 
(P^jta, V £,m ^ 3 in (jvTf)) . consider the effect of 
external magnetic fields on charged quarks produced 
at the early stage of the heavy-ion collisions. More- 
over, in the formalism presented in [2l| . in contrast to 
the dispersion relations presented in (|V.17I) . the en- 
ergy dispersion relation of neutral pion is unaffected 
by the external magnetic field. Using (|V.17[) . and in 
analogy to (|V.15[) . the pions pole masses are defined 

by 



Re[(M|)nTi(M|)i 2 ] 



Re [(J" 00 

Re {Ml) 
Re (-F°°)33 



1/2 



iiTi(^" 00 )i2]J 

-,1/2 

33 



(V.18) 



In particular, the screening mass and the refraction 
index of neutral pions in the i-th direction are given 

by 



,00 



,00 



and 



,(■0 



Re (J 7 '") 



33 



Re (J™) 



33 



1/2 



(V.19) 



respectively. In this paper, we will focus on the In- 
dependence of the mass and refraction index of neutral 
pions at fixed (i and finite eB. The study of the effect 
of constant magnetic fields on charged pion masses 
and refraction indices will be postponed to a future 
publication [30] . 

Let us start with the case T = fi = eB = 0. Using 
the numerical results from (|V.6[) and (|V.7j) . in this 
case, the cr-meson mass and refraction index are given 
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FIG. 6: The diagonal elements of .T MM , fi = 0, • • • ,3 matrices are plotted as functions of T £ [0,400] MeV at vanishing 
chemical potential (jj, = 0) and for eB = 0.03,0.2,0.3 GeV 2 . The identities from ((V773j) for i = 3, T 11 = T 22 and 
J" 00 = — J-" 33 , are explicitly demonstrated in these plots. 



by 

m a ~ 824.3 MeV, and u 
and therefore 

= m a , Vi = 1,2,3. 



a) fi = eB = 



1, (V.20) 



(V.21) 



Similarly, the 7r-meson mass and refraction index at 
T = n = eB = read 



137.7 MeV, 



and 



,00 



V £, i = 1, 2, 3, and therefore 

m« =m w „ W,i=l,2,3. 



1, (V.22) 



(V.23) 



At (T ^ 0,^ = eB = 0), is given by (|VTl5)) . 
Similarly, according to (|V.14|) . the pion masses m# 
are defined by 



Re (T 00 )^ 



1/2 



W= 1,2,3. (V.24) 



Because of the identity (|V.9|) . which is still valid at 
(T 7^ 0, [i = eB = 0), the masses of if = (tti, 7r 2 , 7r 3 ) 
are degenerate, as in T = case [see (|V.22I) ]. In Fig. 
Efa), the T-dependence of to ct and is plotted for 
/i = eB = (black solid line for m a and red dashed 
line for mj). Here, the T-dependence of the coef- 
ficients and T 7 ^ from Fig. [3] is used. We have 
also plotted the T-dependence of the constituent mass 
m = mo + (To in Fig- EI a ) (dotted line) . Comparing 
these curves, it turns out that, as expected, the mass 
degeneracy of a and tt meson masses occurs in the 
crossover region T > 220 MeV. To compare the result 
presented in Fig. [TJa), with the recent results for 
and m a , presented e.g. in (44{, we have set F%" = 
in (|II.10[) . and determined the pole masses of neu- 
tral mesons and the chiral condensate using the same 
method as presented in this paper. The numerical re- 
sults for neutral meson masses and chiral condensate 
for T 7^ 0, fi = eB = and vanishing F^ are plotted 
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FIG. 7: (a) The T-dependence of a and n mesons masses 
is demonstrated at fi = and for vanishing magnetic field 
(solid line for m a and dashed line for m^). Comparing 
these curves with the T-dependence of the constituent 
quark mass m — mo + oq (dotted line), shows that a mass 
degeneracy between m a and mj occurs in the crossover 
region at T > 220 MeV. (b) The T-dependence of the pole 
masses of neutral mesons is plotted for the case when F^" 
in (III.10|l vanishes. The numerical results for m CT and m# 
are in good agreement with the results recently presented 
in 0. 



in Fig. [7|b). As it turns out, only m a changes relative 
to the case where F^" ^ [see Fig. Efb)]. The nu- 
merical results are in good agreement with the results 
presented in [44 |. 

As concerns the screening mass and refraction in- 
dex of if mesons at (T ^ 0, fj, = eB = 0), we use 
the results of Fig. [3j and in analogy to the defini- 
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FIG. 8: The T-dependence of (panel a) and m# (panel 
b) is plotted for eB = and at fj, = 0, 120 MeV. As in the 
case of \i = 0, at \x = 120 MeV, the neutral meson masses 
are still degenerate in the crossover region T > 220 MeV. 



tions (fVU5l) and (fVU6|) define the screening mass and 
refraction index of n mesons by 



(*) 



Re (T l 



where 



1/2 



Re (T 00 ) e 



W,i = 1,2,3.(V.25) 



Using the definitions (|V.16|) and (IV.25|) . and the nu- 
merical results of G 00 as well as (T 00 )^^ from Fig. [3] 
at (T ^ 0, fJ, = eB = 0), the T-dcpendence of the 
screening mass and refraction index of (a, ff) mesons 
can be determined for all directions i = 1,2,3. As it 
turns out, as in T = case, we have 



«« = 1 

(i) 



1. 



and therefore 

(0 
(0 

Tfl^g - 



Mi = 1,2,3, 
W,i = 1,2,3, 



Wi = 1,2,3, 
W,i = 1,2,3, 



(V.26) 



(V.27) 



for the whole interval T £ [0,400]. These results are 
compatible with the identities (jV.101) . The fact that 
Un t = 1 seems to be in contradiction with the results 
from [ll| Hi) , where it is shown that at finite temper- 
ature because of different pion decay constants in the 
spatial and temporal directions, f s and ft, at finite 
temperature, the refraction index u = jr- appearing 



in the energy dispersion relation oj 2 = u 2 (p 2 + in 2 ) 
is smaller than one. Note, however, that in [TH, [24[, 
the pions are self-interacting and f s and ft receive In- 
dependent contributions from one-loop pion-self en- 
ergy diagram, that includes a (tt 2 ) 2 vertex. In con- 
trast, the pions considered in the present paper are 
free. 

Let us also notice that the above results are still 
valid at non- vanishing /i and for vanishing eB. In Fig. 
HI we have compared m a and m ni , I = 1,2, 3 at \i = 
with their values at fj, = 120 MeV for vanishing eB. 
Small deviations from their value at \x = appears 
for m a (panel a). For m n ^,£ = 1,2,3 (panel b) the 
difference between at [i — and [i — 120 MeV 
becomes larger with increasing temperature. As it 
turns out, at non- vanishing chemical potential, the 
degeneracy of the pion masses m 7Tl is still valid for all 
£=1,2, 3. Moreover, m a and toj are also degenerate 
in the crossover region T > 220 MeV for [i = 120 
MeV, as in the fj, = case. 

At finite T and for non- vanishing magnetic fields, 
the pion masses are not degenerate, i.e. we have 
m 3r + ^ m n - 7^ m^o [see (IV.11|) and (|V.12[) and 
the definitions of m 7T ± and m^o from (|V.18j) ]. In 
this paper, we will focus on the T-dependence of m a 
and TTvo. In Figs. |UJa)-|ni[c), the T-dependence of 
(m ff) m T o) masses are plotted for eB = 0.03,0.2,0.3 
GeV 2 and at /i = 0. The expected degeneracy of 
TYLfj and m,o mesons in the crossover region can be 
observed in all the plots of Fig. |H1 However, as it 
turns out, the overlap interval depends on eB for 
fixed fi. Denoting the minimum temperature for 
which the overlap interval starts with T Q , then for 
eB = 0.03 GeV 2 we have T ~ 210 MeV, whereas 
for eB = 0.2,0.3 GeV 2 , T are given by T a ~ 220 
MeV and T c ~ 240 MeV, respectively. 

In Fig. 1101 we have compared the T-dependence of 
the masses of a and ir° mesons for fixed fi = 120 MeV 
and various eB — 0.03, 0.2, 0.3 GeV 2 . As it turns out, 
at temperature below (above) the crossover region, 
the CT-meson masses increase (decrease) with increas- 
ing the magnetic field strength. This qualitative be- 
havior of the T-dependence of m a for various eB ^ 
is comparable with the results presented in [9( (see 
Fig. 3 in j9J). The difference arises from the fact 
that, in contrast to the present paper, the quantum 
fluctuations of a- mesons is considered in Q. And, in 
contrast to the present paper, the contribution of 
appearing in (|II.10|) is not considered in 

Using the definitions of the directional refraction 
index of neutral mesons, iiy and u^d, from (|V.16j) 
and (jV.19l) , as well as the T-dependence of , fx = 
0, • • • , 3 and JW, n = 0, ■ ■ • , 3 from Figs. EJandEl the 
T-dependence of the transverse (i = 1,2) and longi- 
tudinal (i = 3) refraction indices of a and ir° mesons, 
are plotted in Fig. [TTJ From Q 00 = —Q 33 as well as 
(J-oo) 33 = -(J- 33 ) 33 in (TVU21 as well as (TVU31) . the 
refraction index of neutral mesons in the longitudi- 
nal direction turn out to be equal to unity, indepen- 
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a) eB = 0.03 GeV 2 , g = MeV b) eB = 0.2 GeV 2 , g = MeV c) eB = 0.3 GeV 2 , /j = MeV 
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FIG. 9: The T-dependence of m CT and m^o is plotted for eB = 0.03, 0.2, 0.3 GeV 2 and at fj, = MeV. 




FIG. 10: The T-dependence of m„ (panel a) and m^o 
(panel b) is plotted for eB = 0.03, 0.2, 0.3 GeV 2 at a = 120 
MeV. 



dent of T and fj, (see the horizontal red dashed line in 
Fig. [TT]). In contrast, the relations G 11 G 00 as well 
as (^ rll )33 + (^ 00 )33 horn (TVTl) as well as dVT3|) . 
lead to it^ ^ 1 as well as u l ^ 1 for i = 1,2. In 
Fig. [11] the T-dependence of the transverse and lon- 
gitudinal refraction indices of free and neutral mesons 
are plotted for eB = 0.03, 0.2, 0.3 GeV 2 and ^ = 
MeV. As it turns out, in the presence of constant mag- 
netic fields, the transverse refraction indices of neutral 
mesons are always larger than unity. Moreover, the 
transverse refraction index of a (tt ) meson decreases 
(increases) with increasing temperature. Note that 
transverse refraction indices of neutral mesons de- 
crease with increasing the strength of the background 
magnetic fields. It is interesting to add the effect of 
meson fluctuations to the above results and recalcu- 
late the T-dependence of longitudinal and transverse 



FIG. 11: The T-dependence of the transverse and longitu- 
dinal refraction indices of a (panel a) and 7r° mesons (panel 
b) is plotted for various eB. The longitudinal refraction 
index of neutral mesons is equal to unity and independent 
of T (red dashed lines). The transverse refraction index of 
neutral mesons decreases with increasing the strength of 
eB. 



refraction indices of neutral mesons for non- vanishing 
eB and fi. 

Using the definition of the screening masses mS 
from (|VTT5|) and rrS from (|VTT9|) . we arrive at the In- 
dependence of the screening masses of neutral mesons 
at fi = and for fixed eB = 0.03, 0.2, 0.3 MeV. Since 
in the longitudinal direction (i = 3), the directional 
refraction index of a and 7r° mesons is equal to unity, 
the screening masses of the neutral mesons in this di- 
rection are the same as their pole masses m a and m^o . 
In Figs. [F^ancirmi the T-dependence of the screening 
masses of a and ir° mesons in the transverse (i = 1, 2) 
and longitudinal (i = 3) directions with respect to the 
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a) eB = 0.03 GeV 2 , g = MeV b) eB = 0.2 GeV 2 , ^ = MeV c) eB = 0.3 GeV 2 , g = MeV 
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FIG. 12: The T-dependence of the screening mass of a mesons, from (|V.15|I in the transverse (i = 1,2) and 

longitudinal (i — 3) directions is plotted for various eB. As it turns out, the screening mass of the a-meson in the 
longitudinal direction m& ,i = 3 is the same as its pole mass rria (see m CT in Fig. [5J). 
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FIG. 13: The T-dependence of the screening mass of 7v° mesons, m o, from (|V.19|) in the transverse (i — 1,2) and 
longitudinal (i — 3) directions is plotted for various eB. As it turns out, the screening mass of the 7r°-meson in the 
longitudinal direction rn^l,i = 3 is the same as its pole mass m^o (see m n o in Fig. [9]). 



direction of the external magnetic field are demon- 
strated. As it turns out, the screening mass of a and 
7T° mesons in the transverse directions are for all fixed 
eB always smaller than the screening masses in the 
longitudinal direction. In Figs. IT47 a) and (b), we have 
compared the screening masses of a and tt mesons 
in the transverse directions (i = 1,2), respectively. 
Comparing with the plots of Figs. HUT a) and (b), it 

turns out, that, in contrast to m },i = 1,2, the be- 
havior of m&,i = 1,2, by increasing the strength of 
the magnetic field is different from that of m a . And, 
whereas m&\i — 1,2, increases, in general with eB, 
fn^l,i = 1,2, decreases with increasing the strength 
of the background magnetic field. 

At this stage a remark concerning the effects of 
stronger magnetic fields, eB > 0.4 GeV 2 , is in order. 
In Fig. [15l the squared mass of neutral 
is plotted for eB = 0.3, 0.5, 0.7 GeV 2 (or equivalently 
eB ~ 26 ml, 37 m 2 with m v = 138 MeV). As it 
turns out, for eB = 0.5 GeV 2 (eB = 0.7 GeV 2 ), in 
the regime of T < 250 MeV (T < 320 MeV), the 
squared mass of neutral pion is negative. This the pi- 
ons are tachyonic. As we have mentioned before, for 
eB > 0.5 GeV 2 , only lower Landau levels contribute 
to ?7i 2 o [see Footnote 10]. The fact that in the regime 
of LLL dominance and at relatively low temperature 
tachyonic modes appear, is in consistency with the 



recent results presented in [5^]. Here, it is shown, 
that at sufficiently low temperature and in the LLL 
approximation tachyonic instabilities appears in the 
NJL model in 2 + 1 dimensions. The tachyonic in- 
stabilities appearing in m 2 from Fig. [TJ] is another 
example of the appearance of these instabilities at low 
temperature and strong magnetic field in 3 + 1 dimen- 
sional NJL model. 



VI. SUMMARY AND CONCLUSIONS 

In this paper, we studied the effects of uniform mag- 
netic fields on the properties of free neutral mesons, 
a and 7r°, in a hot and dense quark matter. The aim 
was, in particular, to explore possible effects of a back- 
ground (constant) magnetic field on the temperature 
dependence of the pole and screening masses as well as 
the directional refraction indices of these mesons. To 
do this, first, using an appropriate derivative expan- 
sion up to second order, the one-loop effective action 
of a two-flavor NJL model at finite (T, //, eB) includ- 
ing a and ir mesons is determined. Then, using the 
formalism, presented in Sec. HH the masses and re- 
fraction indices of these composite fields are computed 
from their energy dispersion relations. 

As it turns out, the one-loop effective action of this 
model consists of two parts, the effective kinetic part, 
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FIG. 14: The T-dependence of the screening mass of o 
(panel a) and ty° (panel b) mesons in the transverse di- 
rection is plotted for eB = 0.03, 0.2, 0.3 GeV 2 at ^ = 
MeV. 
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FIG. 15: The T-dependence of the squared mass of neutral 
pion is plotted for eB = 0.3, 0.5, 0.7 GeV 2 . 



including non-trivial form factors, and the effective 
potential part, from which we explored in Sec. IV Al 
the complete phase portrait of the model in T — fi, 
T — eB and /i — eB planes for various fixed eB,/i 
and T, respectively. Here, we have mainly reviewed 
the results previously presented in Q for a two-flavor 
NJL model including mesons and diquarks. We have 
shown that the magnetic catalysis of dynamical chi- 
ral symmetry breaking affects the phase portrait of 
this model in two different ways: i) The type of the 
chiral phase transition changes from second to first 



order in the presence of constant magnetic fields, and 
ii) the transition temperatures and chemical poten- 
tials from the chiral symmetry broken to chirally sym- 
metric phase increase, in general, with increasing the 
strength of the external magnetic fields. Only at low 
temperatures T < 50 MeV and high chemical poten- 
tials 280 < /i < 340 MeV and for weak magnetic fields 
eB < 0.2 GeV 2 , the transition temperature decreases 
with increasing the strength of eB. This is related 
to the phenomenon of inverse magnetic catalysis, dis- 
cussed in [8J, [27j . 

In the rest of the paper, we mainly focused on the 
kinetic part of the one-loop effective action. Using the 
formalism originally presented in [2j for a single fla- 
vor NJL model, and generalizing it to a multi-flavor 
system, in Sec. [H] we have determined, in Sec. HVl 
the nontrivial form factors and squared mass matrices 
corresponding to neutral mesons at finite (T, /i, eB), 
up to an integration over ^-momentum and a sum- 
mation over Landau levels. They are then performed 
numerically in Sec. |V1 where, in particular, the In- 
dependence of the form factors and squared mass ma- 
trices of the neutral mesons are presented for several 
fixed magnetic fields and zero chemical potential. Us- 
ing these quantities, we have eventually determined, 
the T-dependence of the pole and screening masses as 
well as the directional refraction index of a and ir° 
mesons for fixed magnetic fields and at vanishing as 
well as finite chemical potential. 

Because of the assumed isospin symmetry, imply- 
ing m u = rrid, charged and neutral meson masses 
are expected to be degenerate for vanishing magnetic 
fields and at zero temperature and chemical poten- 
tial. However, as it turns out, this degeneracy breaks 
down in the presence of constant magnetic fields, so 
that we have m n o ^ ^ m^- , even at zero (T, /i). 
This effect is mainly because of the dimensional reduc- 
tion from D to D — 2 dimensions in the presence of 
constant magnetic fields, which affects the dynamics 
of a fermionic system in the longitudinal and trans- 
verse directions with respect to the direction of the 
external magnetic field. As a consequence, directional 
anisotropy in various quantities corresponding to the 
particles in the presence of a uniform magnetic field 
is implied. 

In the present paper, we have only studied the 
T-dependence of the masses of neutral mesons for 
fixed magnetic fields and chemical potentials. The 
T-dependence of charged meson masses at finite eB 
and /x, will be presented elsewhere [3Cj]. As concerns 
the a- meson mass, m a , the expected mass degeneracy 
with the mass of neutral pions, m^o , in the crossover 
region, T > 220 MeV, is observed for various fixed 
eB and \i. Moreover, as it turns out, m^o decreases 
with increasing the strength of the magnetic field. In 
contrast, m a increases with increasing eB only at low 
temperature T < 220 MeV, while it decreases with 
increasing eB in the crossover region, T > 220 MeV. 
This qualitative behavior is consistent with the re- 
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suit previously presented in [9] in the framework of a 
Polyakov-Quark-Meson model in 3 + 1 dimensions. 

As concerns the refraction indices of neutral 
mesons, it turns out that in the presence of constant 
magnetic fields, the longitudinal and transverse refrac- 
tion indices with respect to the direction of the exter- 
nal magnetic field are different. Moreover, whereas 
the longitudinal refraction index of neutral mesons 
is equal to unity, their transverse refraction index is 
larger than unity. The observed anisotropy in the re- 
fraction indices of neutral mesons is because of the ex- 
plicit breaking of Lorentz symmetry in the presence of 
constant and uniform magnetic fields. The anisotropy 
observed in the directional refraction index of neu- 
tral mesons is also reflected in their screening masses, 
which arc different in the longitudinal and transverse 
directions with respect to the direction of eB. Ac- 
cording to their definitions, and because of the above 
mentioned results for directional refraction indices in 
finite eB, the screening masses of the neutral mesons 
in the longitudinal direction are the same as their pole 
masses, while in the transverse direction, independent 
of T and fx, their screening masses are always smaller 
than their pole masses. They increase with increasing 
temperature at a fixed eB and fx. Moreover, whereas 
the screening mass of a in the transverse direction in- 
creases in general with the strength of the background 
magnetic field, the screening mass of 7r°, in the same 
direction, decreases with eB. 

It is worth to note that the results obtained in this 
paper, showing qualitatively the effect of strong mag- 
netic fields on the properties of neutral mesons in a 
hot and magnetized quark matter, can, apart from the 
physics of magnetars, be also relevant for the physics 
of heavy ion collisions at RHIC and LHC. As it is 
known from [IH, [r| , magnetic fields are supposed to 
be produced in the early stage of non-central heavy- 
ion collisions, and, depending on the initial conditions, 
e.g. the energies of colliding nucleons and the corre- 
sponding impact parameters, they are estimated to 
be in the order eB ~ 1.5 m 2 (eB ~ 0.03 GeV 2 ) at 
RHIC and eB ~ 15 ml (eB ~ 0.3 GeV 2 ) at LHC 
energies. Although the created magnetic field is ex- 
tremely short-living and decays very fast, it can affect 
the properties of charged quarks produced in the ear- 
liest stage of heavy-ion collisions. The way we have 
introduced the magnetic fields in, e.g., (|IV.1[) . where 
the external magnetic field interacts essentially with 
charged quarks, opens the possibility to describe qual- 
itatively the effect of external magnetic fields on neu- 
tral mesons built from these magnetized and charged 
quarks. Note that neutral mesons, by themselves, 
have, because of the lack of electric charge no inter- 
action with the external magnetic fields. Thus, the 
method used in the present paper, is in contrast to 
the method used in [2(J [U| , where the external mag- 
netic field interacts only with charged pions appearing 
in a magnetized chiral perturbative Lagrangian. 

Being motivated by these facts, we mainly focused, 



in this paper, on the effects of weak and intermediate 
magnetic fields, eB < 0.3 GeV 2 . In Fig. [TJJ however, 
we have plotted the squared mass of neutral pion as a 
function of temperature for eB — 0.5, 0.7 GeV 2 . Here, 
we have shown that at low temperature and for strong 
magnetic fields, where LLL approximation is reliable, 
m 2 o becomes negative. The appearance of these kind 
of tachyonic instabilities at low temperature and in 
the presence of strong magnetic fields is recently ob- 
served in 29] in the framework of an NJL model in 
2 + 1 dimensions, which has application in condensed 
matter physics. Our results are consistent with the 
main conclusions presented in [29| . 

Let us also notice that the model used in the 
present paper can be extended in many ways, e.g. by 
improving the method leading to the kinetic coeffi- 
cients and mass matrices of the mesons using func- 
tional renormalization group (RG) method, which is 
recently used in [j| [ljj EH . 
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Appendix: Dimensional Regularization of (|III.13[) 

In this appendix, we will use an appropriate di- 
mensional regularization to regularize the (T, /i)- 
independent part of the effective potential 

n$(m;eB,T = n = 0) 

00 f-\-OG J 

= - 3 E i^i E«w S 1 ^ (A.i) 

9 e{§,-±} p=o J -°° 

appearing in pil.131) . Here, E q is given in (jIII.101) . 
Using the definition of a p = 2 — S p o, we get 

fi$(m; eB,T = fi = 0) = -3 ^ \qeB\ 

x £J £ - E q (p = 0)) ■ (A.2) 
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The above integral can be dimensionally regularized 
using 



+00 id 



d d p 
(2n) d 



+ P 2 Y 
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Setting a = —1/2, e? = 1 — e, with e a small and 
positive number, we arrive first at 

^V(m;eS,T = A1 = 0)^ 3rM + |) 



£ i**i a £ 



4tt 2 
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(aj g + p) 



-1+4 
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where x 9 = 2 i^bi ■ Replacing the sum over the Lan- 
dau levels p with the generalized Riemann-Hurwitz 
(■-function [33], £ (s, a) = £)£i ( a we § et 



n$(m;eB,T = M = 0) = ^ J] (2|geB|) 
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Expanding the above expression in the orders of e up 
to 0(e) and eventually taking the limit e — » 0, we 
arrive at 

n^(m;eB,T = // = 0) 



3e 



(l-7B)(l + 6a# 



- x g lnx g - 2C'(-l,x g ) 
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(A.6) 



Here, we have used the polynomial expansion of 
£(— 1, x g ) = —5 (i — x g + a; 2 ) and the notation 



_ d((s,X q ) 



ds 



(A.7) 



In 1A.6|i . -fE — 0.577 is the Euler-Mascheroni con- 
stant. To eliminate the divergent term, proportional 
to e" 1 in ([AH]) , we use the method introduced in 33], 

and add/subtract to (m; eP, T = (j, = 0) the con- 
tribution of the vacuum pressure 



2\!/2 



(A.8) 



where N c and A/ are the number of colors and flavors, 
respectively. But before doing this, we will first bring 
Pn in an appropriate form. Using (|A.3|) with a = 
— 1/2, setting d = 3 — e, and eventually expanding the 
resulting expression in the orders of e up to 0(e), the 
vacuum pressure, Pq, can be brought in the form 



Pn 



lim 



N c N f m A /(-3 + 2te) 1 km 2 
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Replacing, according to the definition of x q , to 2 with 
rn 2 = 2|qeP|xg, and Nf with a summation over we 



arrive at 
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where N c = 3 is chosen. Equivalently, Pn can be eval- 
uated using a sharp cutoff A [33| , 



P, 
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4 , A + VA 2 + m 2 



to In 



-A(2A 2 + m 2 )v / A 2 + m 2 



(ATI) 



Adding and subtracting P to from (|A.6[) . we 
finally get 



nil\m;eB,T = fj, = 0) 



2tt 



' £ |geB| 2 J C' (-1, x q ) + ^ + ^(1 - lnx, 



+ 4^ ™ 4ln 



A + VA 2 + m 2 



A(2A 2 + m 2 )VA J + m 2 j + a; 9 independent terms. (A. 12) 
I 



Since fi£j (to; eP, T = /1 = 0) is a part of the effec- and we are only interested on the minima of this po- 
tive potential in the gap equation with respect to to, 
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tential, we have neglected the x q (or equivalently m) action, the full effective action of a two-flavor magne- 
independent terms in (|A.12[) . Adding the tree level tized NJL model is given by (|III.14[) . 
and the temperature dependent parts of the effective 
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